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ABSTRACT 

Canonical polar cap cascade models involve curvature radiation (CR) or inverse Compton scattering 
(ICS) of the primary particles and synchrotron radiation (SR) of the higher generation pairs. Here we 
modify such a cascade picture to include the ICS of the higher generation pairs. In such a "full-cascade" 
scenario, not only the perpendicular portion of the energy of the pairs goes to high energy radiation via 
SR, but the parallel portion of the energy of the pairs can also contribute to high energy emission via 
ICS with the soft thermal photons from either the full neutron star surface or the hot polar cap. The 
efficiency of converting particles' kinetic energy to radiation by ICS is very high if the scatterings occur 
in the "resonant" regime. As a result, almost 100% of the energy input from the pulsar inner accelerators 
could be converted to high energy emission. An important output of such a scenario is that the soft tail 
of the ICS spectrum can naturally result in a non-thermal X-ray component which can contribute to the 
luminosities observed by ROSAT and ASCA. 

Here we present an analytic description of such a full polar cap cascade scenario using the recursion 
relationships between adjacent generations following the approach first proposed by Lu et al. (1994), 
but develop it to be able to delineate the complex full-cascade process. The acceleration model we 
adopted is the space-charge-limited flow model proposed by Harding & Muslimov (1998). We present 
the theoretical predictions of the 7-ray luminosities, the thermal and non-thermal X-ray luminosities 
for the known spin-powered X-ray pulsars (8 of them are also 7-ray pulsars), and compare them with 
the observations from CGRO, ROSAT and ASCA. We estimate the non-thermal X-ray luminosity by 
including all the possible ICS branches contributing to a certain energy band, and estimate both the 
full surface and hot polar cap thermal X-ray luminosities by adopting a standard neutron star cooling 
scenario, and by treating self-consistent polar cap heating in the Harding & Muslimov (1998) model, 
respectively. Our results show that the observed different dependences of the high energy luminosities on 
the pulsar spin-down luminosities, i.e., oc (Lgd)^^^ and ~ lO^'^igd, are well reproduced. We found 
that, for normal pulsars, both the hard (ASCA band) and the soft (ROSAT band) X-ray luminosities 
are dominated by the non-thermal X-rays of ICS origin, although for some pulsars, thermal components 
due to either neutron star cooling or polar cap heating can have comparable luminosities so that they 
are detectable. For the millisecond pulsars, our predicted upper limits of the thermal luminosities due to 
polar cap heating are usually higher than the ICS-origin non-thermal components if there are no strong 
multipolar magnetic field components near the neutron star surface, thus the pulsed soft X-rays in the 
ROSAT band from most of the millisecond pulsars might be of thermal origin. 

Subject headings: pulsar: general - radiation mechanism: non-thermal - X-rays: stars 



1. INTRODUCTION 

Among about 1000 spin-powered radio pulsars detected 
so far, 35 are detected as X-ray sources with 11 of them 
being X-ray pulsars (for recent reviews, see Becker & 
Triimper 1997, update version 1999, hereafter BT97; Saito 
1998, hereafter S98; Saito et al. 1997), and at least 8 
are detected also as 7-ray pulsars (for recent reviews, see 
Thompson et al. 1997; Thompson et al. 1999. The eighth 
one, PSR 1046-58, was recently reported by Kaspi et al. 
(1999)). Though detailed data for these spin-powered high 
energy pulsars are still poor, and some of them even have 
a great diversity of emission features, there do exist some 
empirical laws. One prominent feature is the dependence 
of the observed high energy luminosities on the pulsar 
spin-down luminosities: For 7-rays, it is found that oc 
(Lsd)^/^ oc B/p2 (Thompson et al. 1997); while for X- 



rays, the dependence o{L^{ROSAT) - IQ-^L^d oc B^/P^ 
is found in ROSAT (0.1-2.4 keV) data (BT97). A rough 
dependence of L^{ASCA) oc {Lsd f/'^ oc B^/P^ for ASCA 
(0.7-10 keV) data is reported by Saito (1998, hereafter 
S98) and Saito et al. (1997) based on a relatively smaller 
sample than the ROSAT data. However, in view of the 
large uncertainties such as unknown beaming, hydrogen 
column density, distance, assumptions of spectral shape, 
instrument response, calibration uncertainties, and so on, 
the ASCA data are not inconsistent to the ROSAT data. 

Combining the ROSAT and ASCA data, spectral anal- 
ysis of some of the X-ray pulsars are available. Full- 
surface thermal emission components are identified from 
four pulsars, i.e.. Vela (Ogelman, Finley & Zimmermann 
1993), Geminga (Halpern & Wang 1997), PSR 0656 +U 
and PSR 1055-52 (Greiveldinger et al. 1996; Wang et al. 
1998), usually accompanied by a hard non-thermal compo- 
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nent. As for the small-area hot thermal components, the 
evidence is not robust. Geminga (Halpern & Ruderman 
1993), PSR 0656+14 and PSR 1055-52 (Greiveldinger et 
al. 1996) have been reported to show such components, 
but follow-up data analyses show that the hard compo- 
nents are more likely power-law (Halpern & Wang 1997; 
Wang et al. 1998). Two nearby old pulsars, PSR 1929-1-10 
and PSR 0950-1-19 are reported to be dominated by a sin- 
gle small-area hard thermal component (Wang & Halpern 
1997), but such conclusions are also not solid, since PSR 
1929-1-10 could be equally well fitted by a single power law 
(BT97), and the measurements of PSR 0950-1-19 reported 
by Wang & Halpern (1997) are heavily contaminated by 
another source in the ASCA field of view (S98). For the 
millisecond pulsar PSR J0437-4715, two-component spec- 
tral fits using the ROSAT plus EUVE data require a hot 
small-area thermal component (Becker & Triimper 1993; 
Halpern, Martin & Marshall 1996; Zavhn & Pavlov 1998), 
although a simple power-law (BT97) or a broken power- 
law (Becker, Triimper 1999) fit is also acceptable. Among 
others, 7 ROSAT (BT97) and 1 ASCA (S98) pulsars are 
identified to show power-law spectral features, which sug- 
gests a non-thermal origin. The detected X-ray photons 
from the rest of pulsars are too few for spectral fits, but 
they are commonly thought to be of non-thermal origin as 
well, due to the close relationships of their X-ray luminosi- 
ties to the spin-down luminosities (e.g. BT97; S98). 

Two general classes of models have been proposed for 
high-energy pulsars. The polar cap models (Harding 1981; 
Daugherty & Harding 1982; 1994; 1996, hereafter DH96; 
Sturner & Dermer 1994; Sturner, Dermer & Michel 1995, 
hereafter SDM95) interpret pulsar high energy emission 
as curvature radiation (CR) or magnetic inverse Compton 
scattering (ICS) induced pair production cascade in the 
polar cap region; while the outer gap models (e.g. Cheng, 
Ho & Ruderman 1986a, b; Romani & Yadigaroglu 1995; 
Romani 1996; Zhang & Cheng 1997) interpret pulsar high 
energy emission as fan beam radiation originating in the 
large vacuum gaps in the outer magnetosphere beyond the 
null charge surface. In the early polar cap cascade pictures, 
the observed 7-emission is interpreted as the primary CR 
and the synchrotron radiation (SR) from the higher gener- 
ation pairs (Daugherty & Harding 1982, 1994, 1996; Zhao 
et al. 1989; Lu & Shi 1990; Lu, Wei & Song 1994 here- 
after LWS94; Wei, Song & Lu 1997, hereafter WSL97). 
Recently, the importance of the ICS with the soft thermal 
photons from the NS surface by the primary particles is no- 
ticed. It is found that ICS is not only a significant energy 
loss (Kardashev et al. 1984; Xia et al. 1985; Daugherty 
& Harding 1989; Sturner 1995) mechanism, but is also 
a mechanism of 7-ray emission (Sturner & Dermer 1994; 
SDM95) and to ignite pair-production cascades (Zhang & 
Qiao 1996; Zhang et al. 1997a, b). However, Harding & 
Muslimov (1998, hereafter HM98) argued that the asym- 
metry of upward versus downward ICS induced pair cas- 
cades make a ICS-controlled accelerator usually unstable 
near the stellar surface in the space-charge-limited flow ac- 
celeration model, so that a CR-controUed accelerator could 
be stablized at a height of 0.5-1 stellar radii above the 
NS surface when the energy loss efficiency of ICS is less 
than that of CR. As a result, the CR- induced cascades 

^This is true for most of pulsars, but is not the case for the pulsars 
2 and corresponding discussions. 



picture (DH96) could still be maintained, although ICS of 
the primaries will also contribute to the final 7-ray spec- 
trum. The high-altitude accelerator naturally matches the 
extended polar cap cascade model (DH96). 

In such a polar cap model, it is predicted that there is a 
low-energy turnover at the blueshifted local cyclotron fre- 
quency in the high energy spectra of pulsars (Harding & 
Daugherty 1999). Such a feature is also noticed by Rudak 
& Dyks (1999) who argued that the non-thermal X-rays 
observed by ROSAT and ASCA pose a challenge to such a 
model, and it could be a discriminator between the rivaling 
classes of models, since the dependence of ^ 10~^I/sd 
could be interpreted in terms of a thick outer gap model 
(Cheng, Gil & Zhang 1998, hereafter CGZ98; Cheng & 
Zhang 1999, hereafter CZ99). However, there is yet an- 
other important point which has been overlooked in previ- 
ous studies of polar cap cascade scenarios. This is the 
possible contributions of the ICS photons produced by 
higher generation pairs to the pulsar high energy radia- 
tion. Even in the ICS-induced polar cap cascade model 
(SDM95), these contributions were not included, since the 
authors found that these ICS photons are not energetic 
enough to pair-produce again though they admitted that 
most of these scatterings occur in the "resonant" regime 
so that they should be important for spectral formation. 
Here we will improve the polar cap cascade picture by 
including the ICS contributions of the higher generation 
pairs. We will adopt an analytic approach to describe such 
a more complicated "full-cascade" process and examine 
how significant these contributions are. A straightforward 
conclusion is that, the low energy tail of such ICS spec- 
tra can naturally provide a non-thermal X-ray component, 
which could contribute to the ROSAT and ASCA bands 
observations. We will show that, the full cascade model 
combined with the HM98 acceleration model can not only 
present a clear L-y cx (Lsd)^^^ feature, but the ICS-origin 
non-thermal X-ray components with the combination of 
the thermal components due to polar cap heating can also 
reproduce the rough ~ lO^^isd dependence. The pre- 
liminary results of this study is shown in Zhang & Harding 
(1999). 

The organization of this paper is as follows. In Sect. 
2.1-2.4, we describe in detail the analytic approach to de- 
lineate the cascade processes. Such descriptions are inde- 
pendent of the acceleration model and thus could be used 
in any model when the inner accelerator (gap) features are 
specified. Sect. 2.5 reviews the main features of the accel- 
eration model proposed by Harding & Muslimov (HM98), 
which has incorporated the frame-dragging-induced par- 
allel electric fields and both the CR- and ICS- induced 
cascade processes, and therefore is hitherto the most com- 
plete space-charge-limited flow acceleration model. In 
Sect. 3.1-3.4, we explicitly present our model predictions 
of the total high energy luminosities (mainly the 7-ray 
luminosites) , thermal, and non-thermal X-ray luminosi- 
ties from the known pulsars, and compare them with the 
observations from CGRO, ROSAT and ASCA. Our main 
results are summarized in Sect. 4 with some discussions. 



high magnetic fields. See our calculation results presented in Table 
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2. THE MODEL 

2.1. Resonant scattering condition and the "full- cascade" 

picture 

In our work, it is assumed that there exist two thermal 
emission components from neutron stars: one soft compo- 
nent over the whole surface with the temperature Tg, and 
a hot component at the polar cap region with the tem- 
perature Tfi- As reviewed in the previous section, these 
components have actually been detected from some pul- 
sars (though with some uncertainties), and they are com- 
monly expected by various theories. The full surface ra- 
diation is predicted in various neutron star cooling (e.g. 
Van Riper & Lamb 1981; Nomoto & Tsuruta 1987) and 
internal heating (e.g. Alpar et al. 1984; Shibazaki & Lamb 
1989; Schaab et al. 1999) models. In the outer gap heating 
models, it is supposed that the hard X-rays emitted from 
the hot polar cap may be reflected back to the surface by 
a "pair-blanket" due to the magnetic cyclotron resonant 
absorption, so that they might be reemitted from the full 
neutron star surface with softer energies (Halpern & Ru- 
derman 1993; Zhang & Cheng 1997; Wang et al. 1998; 
CGZ98; CZ99). A hotter polar cap is expected by either 
polar cap heating (e.g. Arons 1981; Harding, Ozernoy & 
Usov 1993; Muslimov & Tsygan 1992; HM98) or outer gap 
heating (again Halpern & Ruderman 1993, etc.). The ICS 
of the relativistic particles with these thermal photons is 
the main topic discussed in this paper. Thus our present 
model is solely contingent upon the existence of the ther- 
mal photon fields in the neutron star vicinity, and can not 
be justified if it turns out that there is no such thermal 
emission at all. Note there should be also some ICS pro- 
cesses with the soft non-thermal X-rays if they do exist, 
and such processes will not be discussed in this paper. 

The so called "resonant scattering condition" (Daugh- 
erty & Harding 1989; Dermer 1990) is essential in our 
"full-cascade" picture. This condition ensures that the 
bulk of the soft thermal photons near the neutron star 
surface with peak energy 2.82kT = 0.24T6(keV) could 
be Lorentz-boosted to the particle's rest-frame local field 
cyclotron energy hu/^ = heBe/mc — 11.6i?e,i2(keV), 
where B^. is the magnetic field strength at the emission 
(scattering) point. More specifically, this condition reads 
-yil- (3fi) ■ 2.82kT ;>huj^, or 



48^1 



[(1 - Pf^)n] 



(1) 



where fj, = cos is the cosine of the incident angle of the 
scattering, and 7 is the Lorentz factor of the electron. For 
the case of actual pulsars with two thermal components, 
we adopt different geometries with 



Ms 



(2) 



for the soft full-surface component, and 
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-2U/2 



' e,6) 



(3) 

for the hot polar cap component, where re,6 is the radius 
of the emission point in units of lO^cm (i.e. the stellar 



radius). Note that the two terms on the right hand of (||) 
are for different emission height regimes, the conjuction of 
which is the height above which the horizon of the emis- 
sion point is greater than the polar cap area. Then the 
resonant condition dll) becomes 



7 k, 7ros 
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(4) 

In the resonant ICS regime, the e-folding mean free 
path of the particles with Lorentz factor 7 (Dermer 1990; 
Sturner 1995) 
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is very small due to the tremendous enhancement of the 
scattering cross section near the resonant frequency, so 
that the particles could efficiently convert their kinetic en- 
ergies into radiation before moving a long distance. Thus, 
as long as the condition (^) is satisfied, the Lorentz factor 
of the particles will keep decreasing all the way down to 
the value when the condition (^ fails. Then, most of the 
kinetic energies of the particles are converted to radiation. 

According to Dermer (1990) and Sturner (1995), there is 
a wide regime of the Lorentz factors of the particles which 
satisfies the resonant condition (^). The resonant scatter- 
ing takes place not only for the primary particles within 
the accelerator, but also for the secondary pairs produced 
above the accelerator. Zhang et al. (1997b) have calcu- 
lated the Lorentz factor evolution above the gap for the 
secondary pairs within the framework of the Ruderman- 
Sutherland (1975, hereafter RS75) type vacuum gap accel- 
eration model. They found that the resonant ICS process 
can reduce the energies of the secondary pairs to much 
lower values than the original ones, which may constrain 
pulsar radio emission frequency to below a certain value. 
However, they did not consider the contribution of these 
energies to the high energy radiation of the pulsars. In 
fact, such secondary ICS processes will also produce high 
energy emission which contributes to the final spectra of 
the 7-ray pulsars. The low frequency tail of such spectra 
can contribute to a non-thermal X-ray emission compo- 
nent. As a result, it is important to modify the conven- 
tional CR(or ICS)-SR polar cap model to include this ICS 
effect. 

We then come to a "full-cascade picture" (Fig.l): the 
primary particles within the accelerator emit primary 7- 
rays, which produce secondary pairs via 7— _B process. The 
perpendicular energies of these secondary pairs are rapidly 
converted to radiation via SR, while most of the parallel 
energies will also be converted to radiation via ICS if the 
resonant scattering condition (^) is satisfied. Thus two- 
branch secondary generation 7-rays are produced. These 
secondary 7-rays may undergo the j — B process again 
to produce higher two-branch generation pairs and 7-rays. 
Such processes go on and on until the 7-rays escape from 
the magnetosphere. 
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Fig. 1. — A Schematic illustration of the full cascade process. The 
coordinates of the emission point are {re,9e), and /3e is the angle 
between the radial and tangential directions at the emission point. 
The wavy lines denote the photons. SR-branch photons follow the 
original direction of the pairs, while ICS-branch photons follow the 
tangential direction. I is the attenuation length of each generation 
photon, and Iq is scaled from the original emission point where the 
photon is emitted tangentially. 

In this paper, we will limit ourselves to an analytic 
approach to understand the complex "full cascade" pro- 
cesses. An analytic method (e.g. LWS94; WSL97) is help- 
ful to present an over-all picture for the pulsar population, 
and can, to some extent, predict some features of the nu- 
merical results. 

We have some remarks on the validity of condition (|4|) 
for millisecond pulsars. Since the cyclotron energies in 
millisecond pulsars are much lower than the thermal peak 
energies, very small incident angles (e.g. /x ~ 1) are re- 
quired to keep the thermal photons Lorentz-down-boosted 
to the cyclotron frequency, and we can get an absolute 
upper limit of Lorentz factors of the particles for resonant 
scattering by adopting /i = 1. Noticing 1 — /? ~ 1/27^, 
this turns out 

7 < lOSTg/Sg. (6) 

When this condition fails, as long as condition (|^) is sat- 
isfied, there are still some off-peak thermal photons which 
are resonantly scattered, but the energy conversion efh- 
ciency is much lower due to the smaller "soft" photon num- 
bers. Even if both (^) and (||) are satisfied, the mean free 
path of resonant scattering (g) is still much longer than 
the normal pulsars due to the tremendous decrease of the 
field strength. Thus ICS energy conversion effici enc y is 
low for millisecond pulsars. As shown later in Scct. |3.3| , for 
the millisecond pulsars, usually the non-thermal soft X- 



rays of ICS-origin are dominated by the polar-cap-heating- 
produced thermal X-rays. 

2.2. 7 — _B pair production and photon escaping energy 

The one-photon pair production process or the j — B 
process plays an essential role in the polar cap cascade 
scenario. A detailed description of the process and ex- 
plicit numerical computations were presented in Harding, 
Baring & Gonthier (1997, hereafter HBG97). Here we 
will adopt some approximations within certain reasonable 
regimes and come to some analytic expressions. 

2.2.1. Attenuation mean free path 

We adopt the simple asymptotic formula of Erber 
(1966). Assuming that the one photon attenuation coeffi- 
cient is constant over the distance interested in, we get the 
mean free path of attenuation, which is simply the inverse 
of the attenuation coefficient (e.g. RS75) 



I 



4.4 h Be 

(e^/hc) mc B± 



■ exp 
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3^ 



^6 = 1. 03^2^ exp ( i-^ 



where 



X 



1 E-y B±^ 

2 mc? Br 



(7) 



(8) 



field, B'^ 



m^c^ /eh — 4.414 x lO^'^ is the critical magnetic 
B±/Bc, B± = Ssin^kB, is the energy of 
the 7-photon in the local rest inertia frame, which should 
be corrected by a redshift factor (see eq.||8^) when com- 
pared with the observations, is E^ normalized by mc^, 
and 9]iB is the angle between the photon and the field di- 
rections. Such an asymptotic formula is accurate enough 
when B' = B/B^ < 0.1 (Daugherty & Harding 1983), 
which is the case for most pulsars in the HM98 model 
since the accelerators have been lifted to higher altitudes 
where strength of the local magnetic fields have dropped 
nearly an order of magnitude with respect to the surface 
fields (see Sect.p:5^. But we should bear in mind that the 
asymptotic formula is unreliable at near threshold regime 
(e-ySin^kB ^ 2), and should be modified to a more accu- 
rate asymptotic form (Daugherty & Harding 1983) 



1.035;^ exp 



3X 



(1 + -^) 



(9) 



or even no 



with 5 = O.42(e^/2sin0kB)''-^S'^ooo38^ 
asymptotic formula at all (HBG97). 

2.2.2. Geometric relations 

The importance of general relativistic effects on one- 
photon pair production attenuation in a Schwarzschild 
metric has been examined by Gonthier & Harding (1994, 
hereafter GH94) and HBG97. They found that curved 
spacetime effects may influence the process in three ways, 
namely 1) curvature of the photon trajectories, 2) the red- 
shift of the photon energy, and 3) the change in the dipole 
magnetic field (the shape of the field configuration and 
the enhancement of the field strength). Their conclusion 
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is that the distortion of the dipole in curved spacetime is 
compensated by the curved photon trajectory effect, but 
the increasing of the field strength in curved space-time 
will enhance the 7 — _B absorption and hence, lower the 
escaping photon energy. The relative redshifts of the pho- 
ton energies between different generations are not large 
since the height differences between different generations 
are usually much smaller than one stellar radius. How- 
ever, the redshift effect should be taken into account when 
comparing the calculated frequencies (in the local inertia 
frame) with the observational frequencies (the observers 
are located at infinity) (see eq.|8^). In our analytic study 
below, we will keep the flat spacetime geometry, i.e. as- 
suming strict dipolar field configuration and straight tra- 
jectories of photon propagation. Neither of these assump- 
tions are precise enough, but a combination of the two 
comes to a good approximation (GH94). We only incor- 
porate an enhancement factor of the field strength due 
to curved spacetime (see eq.[^). The redshift effect is 
only taken into account when comparing theoretical re- 
sults with the observations (see eq.[Q). 

In an exact "dipole" field in flat space-time, suppose a 
particle at (re, 0e) emits a photon along the field line, then 
there is a pure geometric relation between the perpendic- 
ular field strength the photon encounters, i.e. B±, and the 
travel distance of the photon I (Hardee 1977) 



and < ^ = ds/0pc < 1 is the ratio of the field line mag- 
netic colatitude at the surface, 9s, to the polar cap angle 

at surface. Ope = (^)^^^ 

For I ^ Tq (this is the case for most of the generations 
except for the last one which will not pair produce again, 
see numerical results of HBG97, thus is a plausible ap- 
proximation to describe the averaged cascade processes), 
we have Va — r^, Ba — B^, and equation (O) reads 



Bi_ ~ -Bf. 



B, 



I 



4 re/9e Pe 



(15) 



where 



4 



Pe ^ - 9.2 X 10^(cm)r,\/fpi/2e-i (16) 

is the curvature radius of field line at the emission point 
(assuming dipolar configuration). Equation ( p^ ) is identi- 
cal with the approximation used in many previous analytic 
studies (e.g. RS75). 

We now incorporate the correction of the field enhance- 
ment effect in curved spacetime. By taking into account 
the general relativistic effect, the magnetic field measured 
in a local rest frame of reference in curved spacetime reads 
(see Wasserman & Shapiro 1983 and GH94, but adopting 
magnetic moment BpR^/2, however.) 



B± =lBpR 



I sin 6e 



X[re{l+Sin^ f3e)+ I COS /3e], 



(10) 



where Pe is the angle between the field direction and 
the radial direction of the emission point, with tan/3e = 
l/2tan6'e (see Fig.l). Note BpR^/2 (rather than BpR^) 
is the dipole magnetic moment of a uniformly magnetized 
sphere of radius R (Usov & Melrose 1995; HBG97), so 
that Bp is identical to the surface magnetic field Bg, and 
is connected with pulsar period P and period derivative P 

by 

Bp = 6.4 X lO^^GiPpy^^, (11) 

which is twice of the value adopted by Manchester & Tay- 
lor (1977) and Michel (1991). 

Near the polar cap region, one has 0e <C 1 , and equation 
(nq) becomes 



B, 



3 2 3 / 



(12) 



where subscript "a" indicates the "absorption point" , so 
that Ta — re + I and Ba = Bp (^^^ ■ Such a B± has a 
maximum when I ~ re/3 (Hardee 1977), with 



or 



Bi.niax ^ 0-OSBp [-] 9e^ 0.08Be9, 



B, 



(13) 



1.2 X 10''GBe,i2P-^^^r'Jet (14) 



Here Bp, Be, Ba are the magnetic fields at the surface, 
the emission point, and the absorption point, respectively, 

9e — S. (^7^)^^^ (according to the dipolar configuration), 



■^curved — -^r.curvc^ ~^ -^0, curve 
3B„i?^cos9 



3B„R-' sin 6 



(1- ^ 

V r 



,1/2, 



^ln(l-^)+l + i^ 



1 In 1 



1 !i 

2 r 



where r„ 



2GM 



'9 - ~^ 

look the modification of the 



(17) 

is the Schwarzschild radius. Let us over- 
dipole" configuration of this 
formula (the effect is actually almost canceled by the ef- 
fect of the curved photon trajectory, GH94), and compare 
it with the magnetic field expression in the flat spacetime 
(note also adopting BpR^/2 as the magnetic torque) 



Bfiat — B^^tiatf + Bgfiiitt 



B„R^ COS 6 



Bp sin 6 
2r^ 



(18) 



We see the curved spacetime effect will increase the 
strength of the magnetic field by a factor of 



q{r) 



[r2 t \ 

I r. curved V J 



B'i 



.curved 



{r)] 



1/2 



flat 



(r)]l/2 



(19) 



With this correction, equations ( p^ , p^|jl^Jl^ ) should be 
modified by multiplying a factor of qljy Note that the 
factor g(r) is equivalent to the function /(??) introduced in 
HM98, which decreases with r, and converges to unity at 
infinity. For a typical neutron star with M = 1.4M0 and 
R = 10km, we have rg/R ~ 0.4, thus the modification fac- 
tor at the surface is q{R) ~ 1.4, and q{r) ~ (1.18 — 1.25) for 
the height of 0.5-1 stellar radii, which is the interested re- 
gion discussed in this paper (see details in Sect. |2.5D. Then 
the maximum perpendicular magnetic field (eq. |14|]) now 
should be modified by multiplying a factor q{re), which is 



B±,„,ax,curved(nPSR) ~ 1.4 X lO^GBe.iaP 



-1/2^1/2 



C (20) 
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for normal pulsars, and 



B 



_L.max,curvci 



d(msPSR) ~ 5.1 X lO^GBe^gP: 



-1/2 1/2 
3 'e,6 



e (21) 



for millisecond pulsars. Here — Bp (^-^^ is still the 

field strength at the emission point in flat spacetime. 

There is yet another effect, namely the "non-zero pitch 
angle" effect, that will reduce the attenuation length of 
the SR photons, which has been overlooked in some of the 
previous studies. The expressions of B± (eqs. |To| , p^ , [T5t ) 
are derived assuming that photons are emitted along the 
tangent direction of the field line, which is not the case 
for higher generation SR photons which follow the direc- 
tion of the pairs which have a non-zero pitch angle with 
respect to the field lines. Thus equations ( |l0| , [12|jl^ ) have 
in fact underestimated the strength of the perpendicular 
field. However, the geometric relation in ( ]lO| ) could still 
be used by replacing the emission point (re, 9e) with the 
point {re,o, ^e,o) where the 7-rays are emitted tangentially. 
Defining Iq as the distance between (re.o, ^e,o) E^nd [ra, 
9a), and / as the mean free path of the SR 7-rays, i.e. the 
distance between (re, 9e) and (ra, 9a) (see Fig.l for illus- 
tration), again adopting 0e <C 1, the perpendicular field is 
then 



B± — Bat 



which is simply 



'2ra{re+lo-l){'2re + lo-l)' 



B, 



B, 



Pe 



(22) 



(23) 



when I < Iq r, 
ing I to Iq 



This is just equation (|T^) by chang- 
We see the higher the SR generation, the 
stronger the perpendicular fields, since Iq increases with 
generations. This effect will raise th e valu e of x for pair- 
production with the same Iq (see Sect, 2.2.^ and Fig. 2) and 
hence, tend to increase the number of the SR generations. 
Such a "non-zero pitch angle" effect is non-relevant to the 
ICS branch, however, since the ICS photons are emitted 
by the pairs after emitting their perpendicular energies via 
SR, and hence, follow the tangential direction of the field 
line at the emission point. 

2.2.3. How large is x for pair production? 

Equation (|^) can give a self-consistent Iq and x for one- 
photon pair production. It is essential that the parameter 
X does not sensitively depend upon various parameters, 
since small changes in x correspond to exponentially large 
changes of Iq. Such a feature makes some simple analytic 
formulae available in analytic studies (e.g. RS75; LWS94; 
WSL97). In these models, x is adopted as 1/15. However, 
in order to treat the process more precisely, it is still worth 
examining the possible variation of x for different pulsar 
parameters. A smaller Iq can increase x and hence increase 
the total number of generations (i.e. the gene rati on o rder 
parameter, LWS94, WSL97, also see Sect!|j and|U). In 
fact, Zhang et al. (1997b) has adopted x = 0.09,0.12 to 
estimate the parameters of different polar gap modes. 

Using (^, and (|l^), one can get a more accurate 
value of X when i?e, P, r^fi (these three parameters are 
specified for a certain pulsar and an acceleration model). 



Iq (this parameter depends on the 7-ray energy) and ^ 
(depends on the field line) are given 



X--ln-i[0.01SeP-'/'/6' 



-1/2 

e.6 



(24) 



Figure 2 plots the dependence of 1/x on Iq (correspond- 
ing to different E^) for two sets of parameters: one set is 
for a typical normal pulsar with B^ — lO^^G, P = 0.1s, 
re, 6 = 1-8 and ^ = 1; another is for a typical millisecond 
pulsar with Be = 5 x lO^G, P = 0.005s, re,6 = 1-2 and 
^ = 1. Note that different ado ptions of re,6 are based on 
the HM98 model (see Sect.|U). We found that, the pa- 
rameter X has deviated from 1/15 for most of the cases. 
For most of the generations, we have <C 1 (see numer- 
ical results in HBG97), thus x is larger than 1/15. For 
millisecond pulsars, the values of x are even larger. 



X 



1 , 1 


















1 





Log(g 

Fig. 2. — The 1% dependence of the critical value of x for 7 — -B 
pair production (1/x is plotted). The upper two lines are for a typ- 
ical normal pulsar with P = 0.1s, Be, 12 = 1-0, re,6 = 1.8 (thus 
Bp, 12 = 5.8), and J = 1; the lower two lines are for a typical mil- 
lisecond pulsar with P = 0.005s, Be, 12 = 5.0 X 10~*, re, 6 = 1-2 
and 5 = 1 (thus log Bp = 8.9). The dashed lines are for the case of 
'tangentially' emitting, and the solid lines are for the case with an 
estimate correction of the non-zero pitch angle effect. The dotted 
line marks the conventional value x = 1/15. 

The critical x for pair production could be further en- 
hanced by including "non-zero pitch angle" effect as dis- 
cussed above. With (^ rather than (p^, the equation 
d24) now becomes 



-1/2 



X= ^\n-'[OmB,p-^^^lQlo,Qr^^ 



(25) 



At higher generations, Zo.e is larger than Iq. We have 
adopted Iq^q ~ ZIq for normal pulsars and Iq,q = 1.5^6 
for millisecond pulsars to cover this effect. The different 
adoptions reflect different generation numbers in normal 
and millisecond pulsars. The results are also shown in 
Fig.2. 

With the results in Fig.2, we come to some "recom- 
mended" values for the critica l y. W hen calculating pho- 
ton escaping energy (see Sect. 2.2.4), it is appropriate to 
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adopt a smaller x since /g of the last generation 7-rays 
is usually around unity (recall i?max is achieved when 
/ ~ re/3, also see numerical results of HBG97). However, 
when describing the cascade process, a larger x should be 
adopted, since the averaged is much less than unity. In 
our following analytic formulae, we'll adopt 

Xcsc(nPSR) - 1/16, Xcas(nPSR) - 1/12 (26) 

for normal pulsars, and 

Xcsc(msPSR) - 1/12, Xca.(msPSR) - 1/8 (27) 

for millisecond pulsars, according to the results of Fig. 2. 



2.2.4. Photon escaping energy 

Photon escaping energy is the energy below which the 
magnetosphere is transparent. A numerical result of the 
escaping energy of photons is shown in Fig. 2 of HBG97. 
With (H), adopting B±^ as the i3j_,inax described in (20) 
and ( [21^ in which the correction of curved spacetime ef- 
fect has been included, also choosing the recommended 
Xosc in equations ( p6| , p7|) , we get 



£;^,csc(nPSR) ~ 2.0GeVS, 
for normal pulsars and 



-1 pi/2. 



e,6 ? 



Xi 



/16 



(28) 



£;^,esc(msPSR) ^ 73GeVB-ip!f r^ j/'r'Xvi2 (29) 

for millisecond pulsars. Note these relations are quite 
different from equation (36) of Hardee (1977), who has 
adopted the incorrect form of 7-ray absorption formula. 
As pointed out by Zheng, Zhang & Qiao (1998), although 
the general formula of 7-ray attenuation including both 
the perpendicular magnetic fields and the electric fields 
derived by Daugherty & Lerche (1975) is correct, when ap- 
plying it to the situation of pulsars, the contribution of the 
induced electric fields are actually of no importance, since 
in the static frame (rather than the co-rotating frame), the 
emission direction is not strictly along the direction of the 
field line due to Lorentz transformation. Hardee (1977) 
has used the Daugherty & Lerche's formula by assuming 
that the photons are emitted tangently along the field line 
in the static frame, and thus gave a incorrect larger ab- 
sorption coefficient^. WSL97 obtained a similar equation 
(their eq.[10]) as our (|2^), but they did not include the 
general relativistic correction and the non-zero pitch an- 
gle effect. 

2.3. Recursion relations between different generations 

There exist some simple recursion relations between the 
characteristic emission frequency of both the two-branch 
filial-generation (SR and ICS) photons and the parent- 
generation photons. Consider a certain parent generation 
7-ray with dimensionless energy Ci (the origin of this pho- 
ton is not important), for SR branch, the typical energy 
of the filial-generation 7-rays is 



= -7,^+1 sin 6lkB = 7Xei 



where 7i+i — 6^/2 is the Lorentz factor of the filial- 
generation pairs following the direction of the parent pho- 
ton, which impacts the field with an angle of 0kB ■ Define 
K = e^+i/ci as the energy reduction factor of the adjacent 
generations, we have (LWS94; WSL97) 



(31) 



Recalling the recommended Xc 
(|2^j27|), we have 



values in equations 



KSR(nPSR) - 1/16 



for normal pulsars and 

KSR(msPSR) 



1/11 



(32) 



(33) 



(30) 



for millisecond pulsars. Note in (30) we have adopted 
the classical formula to describe the SR process. In the 
strong magnetic fields, however, the SR process should 
be strictly described quantum-mechanically as the transi- 
tions between Landau levels with adjacent energy interval 
of Ai^L ~ TiWg = 11.2i3e.i2keV (for a review, see Hard- 
ing 1991). The comparisons of the quantum and classical 
treatments of the SR process have been presented in Hard- 
ing & Preece (1987) and Baring (1988). The classical de- 
scription is not a bad approximation as long as the Landau 
state number of the pairs is large (e.g. > 20), which means 
that critical SR photon energy is much larger than AE^, 
and the field strength is not too high, i.e., B' < 0.1 (see 
numerical results in Harding & Preece 1987, note this is 
the same condition for the asymptotic formula (|^) to hold) . 
Both conditions are fulfilled in most cases discussed in this 
paper (see Sectj2.5|). 

For the ICS branch filial-generation, the typical energy 
of the emitting (scattering) photon {{i + l)-th generation) 
is determined by the Lorentz factor of the pairs {{i + 1)- 
th generation) after they have emitted their perpendicular 
energy via SR. Suppose such Lorentz factor is 7^+1, y, the 
typical ICS photon energy is simply (e.g. Zhang et al. 
1997a) 

e,+i,ics = '2l^+l,\\Bi (34) 

for resonant scatterings. Note essentially, this typical en- 
ergy depends on the local field strength, but does not de- 
pend on the incident angle of the scatterings. This feature 
can be understood as follows: a soft photon with frequency 
eo which impacts by an incident angle of 6 — cos^^ /i with 
an electron with energy 7|| , will be inverse Compton scat- 
tered to a maximum frequency of = 27j|(l — /3||/i)eo. 
For a resonant scattering occuring, the incident energy 
eo should be Lorentz-boosted to exactly the cyclotron 
energy €3 — B' in the electron rest frame, i.e., eo = 
^'/[7||(1 ~ Thus the angular factors in both ex- 

pressions are exactly canceled in the final result. In other 
words, the electron can "pick up" the right energies of 
the photons from different directions and resonant scatter 
them to contribute to the same characteristic energy (|^ . 



^Note that in Zheng 
(1981), and Daugherty 
directions" in the static frame. 



2 

Zhang & Qiao (1998), there are some mis-comments on the work of Harding, Tademaru, 
Harding (1982), who have correctly treated the absorption problem by including the 
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"aberration of the emission 
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Another point is that in (p^, we have again overlooked 
the possible variation of B^. for different generations. (We 
did not specify Be,i for a certain generation i.) We will 
adopt this assumption throughout the paper, since it is 
sound only except for the very last generation or when res- 
onant scattering condition (Q) fails so that the ICS mean 
free path (^) becomes long. Even so, the correction to our 
results is still tiny. 

The parallel Lorentz factor of the {i + l)-th generation 
pairs is (see Appendix A, also WSL97) 



or 



7^+l.|| 



7»+i 



[l + (7^i-l)sin2M^/^' 



(35) 



Noting that ji^i — 6^/2 (each of the (i + l)-th pair gain 
one half of the i-th generation photon energy), as long as 
7^+1 » 1 (this is the case for most generations), we have 
(7f+i - 1) sin" 0kB ^ (ez/2 sin^kB)' = {x/B'J^, which only 
depends on when we adopt a constant x- Define 



m 



2M 



(36) 



as the portion of the particle energy which goes to the ICS 
branch, from (p5f) we obtain 



1 



[(x/sO^ + i]V2- 



(37) 



This is independent on the generation order i if we ne- 
glect the very small variation of Be in different genera- 
tions. Note this is an essential point, since the energy 
distribution between the two-branch (SR and ICS) filial 
generations is settled for a giv en pul sar in different gener- 
ations of the cascades. Using (pap^, we get 



77||(nPSR) 



for normal pulsars, and 



1 



{13.5B^J, + iy/^ 



77||(msPSR) ~ 1.8 X lO^^Be,9 



(38) 



(39) 



for millisecond pulsars. We see only a small portion of en- 
ergy is left in the ICS branch for millisecond pulsars, since 
Be is much lower. This is the third factor to reduce the im- 
portance of the ICS-branches for the millisecond pulsars. 
Recall the other two factors are: longer mean free path 
(eq.[||) and possible failure of the second resonant condi- 
tion (eq. §]). These are the reasons why the non-thermal 
X-rays are usually not important compared with the ther- 
mal emi ssion for millisecond pulsars, as will be shown in 
Sectjsl. 



We can also define the energy portion going to the SR 
branch. This is just 



VI- 



1 



(40) 



and was denoted as "energy conversion efficiency" in 
WSL97. Similarly, we can derive different 77j^(nPSR) and 
?7j^(msPSR) using ( |3^ ) and (39), respectively. 

With (||), noting 7i-i-i,|| = '7||7i-i-i> and again 7^+1 = 
ei/2, we have 



KICS 



ei+i,ics 



^V\\B'e, 



(41) 



Kics(nPSR) 



0.0227Se.i2 



1 



(13.55-22 + 1)1/2 - 160 ^"^-i" (42) 

for normal pulsars (assuming i?e,i2 ^ 1) and 

Kics(msPSR) ~ 4.1 X IQ-^Blg (43) 

for millisecond pulsars. This means that, except for pul- 
sars with high magnetic fields, ICS branches usually have 
a more significant decrease of the typical photon energy 
than SR branches, especially for millisecond pulsars. 

2.4. Generation order parameters 

The so-called "generation order parameter" was first in- 
troduced by Lu, Wei & Song (LWS94) and improved by 
Wei, Song & Lu (WSL97). In their work, the authors in- 
troduced the non-integer generation order parameter ( to 
describe the conventional curvature-synchrotron cascades 
based on the RS75 acceleration model, and showed that 
the observed 7-ray pulsars tend to have relative large Cs. 
Such a concept can delineate the sketch of the cascade pro- 
cess, and, to some extent, compensate the shortcoming of 
the analytic "recursion" descriptions, i.e., the ignorance of 
the detailed spectral features of different generations. 

In our full-cascade picture, the complicated two-branch 
filial-generation feature makes it impossible to describe the 
whole process with one single generation order parameter. 
But we can still borrow their definitions in our studies. 

Suppose the primary generation 7-rays have a typical 
energy Eq (this is model-dependent, we'll give the expres- 
sions in the HM98 model (eqs.||,|l[) in Sect.^l), there 
will be altogether 



Csr 



log(5esc/-Eo) 
log(KSR) 



(44) 



SR generations before the photon energies are reduced to 
Ecsc- Note the term (-1-1) in ( ^ ) is to ensure C = 1 for 
the primary 7-rays (with energy Eq). is just the gen- 
eration order parameter defined by LWS94 and WSL97. 
In our full cascade scenario, this parameter only describe 
the typical cascade branches with pure SR generations. 
Similarly, we can define another parameter 



\ogiE,sc/Eo) , ^ 

yes - — j — -f ^ H i 

log(Kics) 



(45) 



to describe the pure ICS cascade branches. However, 
neither ( ^4| ) nor (45) can describe the "crossed" gener- 
ations (e.g. ICS branch from the SR parent generation 
or vice versa). We will define some more general genera- 
tion order parameters (e.g. [Q,|^) to construct an an- 
alytic estimate formula of pulsar X-ray luminosities (see 
eqs.[^|4j). For most pulsars, the ICS branch photons 
usually do not pair-produce any more (see Sect. |2.5| , also 
SDM95). Thi s ca n further simplify the description (see 
details in Sect. 3.3). 



ZHANG & HARDING 



9 



2.5. HM98 acceleration model 

Here before, we have not specified an acceleration model. 
The descriptions above are valid for any polar cap accel- 
eration model once Eq is specified. There are generally 
two subclasses of polar cap acceleration models. One sub- 
class is the vacuum gap model first proposed by RS75, and 
modified by Usov & Mehose (1995; 1996) and Zhang et al. 
(Zhang & Qiao 1996; Zhang et al. 1997a,b). Another sub- 
type model is the space-charge-limited flow model (Arons 
& Scharlcmann 1979; Arons 1983), which assumes free flow 
of particles from the surface. Such a model is improved 
by Muslimov & Tsygan (1992) and Muslimov & Harding 
(1997) by incorporating the frame dragging parallel elec- 
tric fields, and finally HM98 imposed the zero E\\ boundary 
condition at the pair formation front (PFF) to treat the 
accelerator more realistically. In reality, the environment 
of different pulsars could vary much from one to another, 
so that the above-mentioned two different types of inner 
accelerators may anchor in different pulsars. For exam- 
ple, there is growing evidence showing that a handful of 
the "drifting" radio pulsars seem to have the RS-type vac- 
cum accelerators in their polar cap region (Deshpande & 
Rankin 1999; Vivekanand & Joshi 1999). However, theo- 
retical arguments, i.e., the so-called binding energy prob- 
lem encountered by the RS75 model (Jones 1985, 1986; 
Neuhauser et al. 1986, 1987), support the space-charge- 
limited flow acceleration scheme. In this paper, we will 
adopt the HM98 acceleration model uniformly for all the 
pulsars. 

No simple analytic description for all pulsars is available 
in the HM98 model due both to the complicated form of E\\ 
within the gap and to the screening details. Nevertheless, 
a general accelerating picture is available, i.e., grows 
approximately linearly if the gap length Sc is smaller than 
the "effective" polar cap radius rp^E (assuming an effec- 
tive star surface Re ^ (1-5 — 2)i?), or will be saturated to 
an almost constant value above Tpc.E if Sc > Tp^.E- Here- 
after we will define the cases with or without field satura- 
tion as the regime I and II, respectively. In principle, both 
CR and ICS of the primary particles should be taken into 
account. However, for most pulsars (except for high B pul- 
sars), as pointed out by HM98, as long as ICS dominates 
in the accelerator, the accelerator itself is unstable, since 
the upward ICS occurs in the resonant regime, while the 
downward ICS occurs in the Klein-Nishina regime. Thus 
the final stable accelerator should be controlled by CR. In 
fact, the "effective surface" i?^; is determined by the cri- 
terion that beyond this height, the ICS energy loss rate is 
less than the CR energy loss rate (see HM98). In regime I, 
the CR-controUed gap length (S'c), the maximum typical 
Lorentz factor of the primary particles (70), as well as the 
characteristic CR photon energy emitted by the primary 

particles, i.e., Eq = could be approximately ob- 

tained using (A3) of HM98, which read (see also eqs. [77,79] 
in HM98, in which there are some minor computing errors) 

Sc{l) ^ 4.2 X lO^cmSp^ff p4/7^i6/7(^Qg^)-3/7^ (4g) 
7o(I) ^ 8.2 X IQ-'B^li'P^/'R'ilicosaY'^ (47) 



and 



respectively, where a is the inclination angle of the neu- 
tron star. For the regime II, we adopt the saturated E\\ 
presented in (A5) of HM98, and obtain the different ex- 
pressions which read 



Sc{ll) ^ 3.2 X 10^cmBj;;j2^^^*i?^2g(cosa)' 



-3/4 



70(11) 



1.4x 10^p-l/4i^|6(cosa)l/^ 



and 



Eo(ll) ~ 0.88(GeV)p-5/4i?ef (cosa)3/4 



(49) 



(50) 



(51) 



The emission height discussed in the previous sections is 
then 



^efi — Res + Scfi 



(52) 



It is worth noting that, the primary particles actually lose 
their energies within and above the accelerator via both 
CR and ICS, and the ICS process has a different typical 
Eq, so that it will also contribute to a component in the 
final 7-ray spectrum and change the generation structure 
to some extent. But such modifications are not prominant. 
Furthermore, the energetics of both the CR and ICS pho- 
tons are from the kinetic energy of the primary particles 
(|4^j5^) , thus the variation of the detailed primary photon 
spectrum does not change much the luminosity predictions 
presented in Sect. 3, where we adopt the primary luminos- 
ity as the total polar cap particle luminosity (see Q). 

Since regime I is limited by regime II, in principle, we 
should adopt the minimum of 70(1) and 70 (II) as the real 
7o (and hence Eq) for a certain pulsar. The critical con- 
dition to separate the two regimes could be obtained by 
equating eq.(p7|) with ([5C|). Thus a pulsar should be in 
regime I if 



</,>-"/^«<//(cosa)^/^«>6.0. 



(53) 



i?o(I) ~ 178(GeV)Pp;?f p-l/l44'/6"(cosa)^/^ (48) 



N3/7 



and in regime II otherwise. Note that this is just a 
rough analytic treatment. Numerically, there should be 
no such distinct boundary, and the two regimes should 
meet smoothly. 

With (^ ) or (51), one can get explicit expressions of 
the generation parameters from (^^ and (^). This com- 
pletes the analytic description of the full cascade process. 
It is worth noting that 1 < Cics < 2 for most pulsars (see 
Table 2). Such an issue was also noted by SDM95, who 
thus suggested that ICS branches will not further produce 
filial-generation pairs. However, for the high field pulsars, 
this is not necessarily the case (also see Table 2). In our 
analytic formula below, without losing generality, we will 
also cover the possible contributions from the ICS-photon- 
produced pairs (see eq.[[77|). 

HM98 has discussed the positron backflow feedback 
quantitatively. Though the complication of the E\\ screen- 
ing process prevents an accurate determination of the frac- 
tion of backflowing positrons, an upper limit of such a 
flow is readily available from the model, which is useful 
to estimate the maximum thermal X-ray luminosities due 
to polar cap heating. The expressions of such maximum 
thermal X-ray lu mino sities are shown in equations (69) 
and (^ (see Sect.^. 
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3. LUMINOSITY PREDICTIONS 

3.1. Gamma-ray luminosity 

Harding (1981) pointed out that by assuming a roughly 
constant potential in the inner accelerators for different 
pulsars (so that the primary particles gain almost the same 
energy 70 for different pulsars) , the total polar cap "lumi- 
nosity" 

Lpc = jomc^Np (54) 
is then proportional to Bp/P^, since the particle flow rate 



A natural implication of the Lpc oc (Lgd)^^^ dependence 
is that the 7-ray emission efficiency, Lpc/isd, will increase 
for the pulsars with lower Lsd- Physically, such an effi- 
ciency can not be greater than unity. Since as Lpc ap- 
proaches Lsd the accelerator should be in regime II (satu- 
rated), we can use (|60|) and (^8|) to get the constraint 



1.4 X 10-^'' R^Bp^uP-^ [cos a) (55) 



itself is proportional to Bp/P^, where n^j ~ |^ is the 
Goldrcich- Julian (1969) number density. This will give 
a natural interpretation of the observed relationship of 
oc (Lsd)^^^ if the polar cap "luminosity" is completely 
converted to 7-ray emission. However, in the canonical 
CR-SR (DH96) or ICS-SR (SDM95) cascade picture, the 
expected 7-ray luminosity Lj is not fully identical to Lpc, 
since only the perpendicular portion of the energies of the 
pairs goes to the next generation radiation. More specifi- 
cally, one has (WSL97) 



L^ (canonical) 



r „CsR 



(56) 



which will deviate from Lpc to some extent. The deviation 
is more significant for the pulsars with small 77_l and large 
CsR (e.g. PSR 1509-58, see Fig.3, Table 1). In the fuU 
two-branch cascade picture described above, it is safe to 
regard 

L^(full) ~ Lpc, (57) 

since the ICS branches "pick up" the "lost" parallel kinetic 
energies of the particles also to radiation, so that almost 
100% of the polar cap "luminosity" is converted to high 
energy radiation. Although not all these energies are con- 
verted to 7-rays (in fact, the non-thermal X-r ay emission 
is also part of this energy budget, see Sect. 3.3 below), the 
majority of this energy output is in the 7-ray band. 

In the HM98 model, the typical energy of the primary 
particles 70 turns out to be very weakly dependent on the 
pulsar parameters (eqs. |47|j5C| ]), especially for the young 
pulsars in regime I (only Geminga and PSR 1055-52 are 
in regime II) . This is in agreement with the assumption of 
Harding (1981). With eqs.(|^,||j55|) and (|7|) or (|^, and 



9.68 X 103°s2_i2 



(58) 



(/ = 10"*^ is the moment of inertia), the 7-ray luminos- 
ity is finally 
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L^(full,I) ~ 9.4 X lO^^erg • a'^ Bl'J^p-^'^''' R%l{cosa) 

= 3.0 X lQ^^B-^\gP^/'R%l{cosaf\UdY/\ (59) 
for regime I, and 

LT,(full,II) ~ 1.6 X lO^^erg • s'^ Bp,i2P-^''^R%,Q{coiiaf'^ 



- 0.5 X 10l6p-l/4^|^^(cosa)5/4(L,d)l/^ 



(60) 



for regime II. This almost reproduces the L-^ cx (Lgd) 
feature (also see Fig.3). 



1/2 



Bp,i2P~'^'^RE%{cosa)-^'^ > 1.65. 



(61) 



Such a line is also close to the pulsar pair formation 
"deathline" (for a detailed discussion of the deathlines, 
see Zhang, Harding & Muslimov (1999)). The luminosity 
formulae (^,^^ are no longer valid when (61) fails. None 
of the observed spin-powered X-ray pulsars studied in this 
paper are beyond this line. 

In Table 1 and Fig.3, we compare the broad-band high 
energy luminosities of the 8 7-ray pulsars with the the- 
oretical predictions both in the full cascade model and 
the canonical CR-SR cascade model. Though the oc 
(Lsd)^^^ feature is reported using 7-ray luminosities above 
lOOKeV (Thompson et al. 1997), here we adopt the lat- 
est luminosity data above leV for the 7 previously known 
7-ray pulsars including the new data of PSR 1055-52 
(Thompson et al. 1999), which also clearly shows the same 
trend, and indicates that most emission is spread in the 7- 
ray band. For the newly discovered 7-ray pulsar PSR 1046- 
58, we adopt the luminosity data above 400MeV (Kaspi et 
al. 1999). As shown in Fig.3, all these data just meet our 
model prediction L-y(full) (cqs.[^,^) well, which is essen- 
tially the polar cap luminosity Lpc (eq.Q). The model 
parameters we adopt are: a = 30° and i?_E,6 = 1-8. It is 
clear that the full cascade model shows a much better fit 
to the observations than the canonical model which pre- 
dicts lower 7-ray luminosities (eq.|56|) (see Fig.3), espe- 
cially for PSR 1509-58. The small differences between the 
full-cascade predictions and the observations are within 
the range of different a and slightly different Re.g- Thus 
our model is in good agreement with the observations. 
Note that the present description is incomplete for the 
high B pulsar PSR 1509-58 (with log Bp = 13.49). In 
such strong fields, the resonant frequency is much higher, 
and the downward ICS of the positrons also occurs in the 
resonant regime (HM98). The anisotropy of ICS processes 
then dissappears and the accelerator should still be located 
at the neutron star surface, with ICS as the controlled 
mechanism. However, photon splitting will also play an 
important role in high B pulsars (HBG97), so that the 
cascade is more complicated and includes both pair pro- 
duction and photon splitting. All these effects may alter 
the luminosity predictions of high B pulsars significantly, 
both for 7-rays and also for X-rays as well. 
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Fig. 3. — Comparison of the observed broad-band luminosities 
for the 7 7-ray pulsars and 7-ray luminosity of PSR 1046-58 above 
400MeV with the model predictions in both the full cascade and 
the canonical CR(ICS)-SR cascade models. Note how the full cas- 
cade model 'picks up' the luminosities which the canonical model has 
'lost', and clearly reproduces the Lj cx (Lgd)"'^''^ feature. 

It is also of interest to examine 7-ray luminosities of the 
millisecond pulsars. According to (pOh, usually the mil- 
lisecond pulsars also have considerable 7-ray luminosities 
(oc (Lgd)^^^)- However, due to their weak field strengths, 
the escaping photon energy is quite high (eq. ) . As a re- 
sult, 7-ray spectra of the millisecond pulsars are expected 
to be very hard, with the lowest characteristic SR energy 
(eq.[f4|) about ISOGeV, and the detectable photon num- 
ber isiewer. Future missions, such as GLAST, may detect 
emission from these sources. 

3.2. Thermal X-ray luminosity 



As reviewed in Sect. 2.1, there are many possible mecha- 
nisms to account for both the full surface and the hot po- 
lar cap thermal emission components from neutron stars. 
Here we will treat the full surface emission by adopting a 
simple "standard" cooling model, and treat the hot spot 
emission using the polar cap heating scenario predicted by 
the HM98 model. 

3.2.1. Full surface thermal emission: neutron star cooling 

Neutron star cooling is the results of many different 
mechanisms, the physical details of some of which are still 
poorly known. Thus there are many uncertainties in differ- 
ent models (see Shapiro & Teukolsky (1983) for a review). 
In our calculation, we use a rough relation 



logT,.6 
logT,,6 



-0.1 log r 
-0.5 logr 



0.37 
2.45 



logr < 5.2 
logr > 5.2 



(62) 



which is obtained by a rough fit to the latest detailed nu- 
merical model of Schaab et al. (1999, their Fig. 3), and is 
consistent with the standard "modified URCA" neutrino 



cooling plus photon cooling model reviewed in Shapiro & 
Teukolsky (1983). Here, r = P/2P is the spin-down age 
of the pulsars. The reasons for us to adopt such a sim- 
ple treatment are as follows: Firstly, such simple cooling 
models are not inconsistent with the observations (BT97; 
Schaab et al. 1999). Secondly, as long as the resonant 
ICS condition (^) is satisfied, our non-thermal luminos- 
ity predictions are insensative to the temperature of the 
soft photons. Thus some extent of inaccuracy in the full- 
surface temperature adoption has little influence on the 
final results. Thirdly, in the polar cap cascade scenarios, 
it is unlikely to form a nearly static pair "blanket" near 
the surface, as assumed in the outer gap models. Hence, 
there should be no "reflected" full-surface thermal emis- 
sion component. Furthermore, the outer gap models did 
not compare their soft thermal emission prediction with 
the observations, so that the assumption can not be justi- 
fied. 

3.2.2. Polar cap heating 



As discussed in Sect.^J], the complication of the 
screening process makes it difficult to treat the positron 
backflow accurately, so that an explicit analytic descrip- 
tion of the polar cap heating is not available. In principle, 



the maximum backflow positron rate 



could be es- 



timated by multiplying the polar cap current Np (eq. [£5[ ) 
with a factor 

(V-E)||/(47r) 



/ 



2Pc 



that^ 



pcj 



(63) 



(64) 



since determines the charge density required for sceen- 
ing, and the factor (1/2) accounts for the readjustment of 
the primary current due to the backflow current. How- 
ever, due to the mathematical complications of the space- 
charge-limited flow model with frame-dragging included 
(Muslimov & Tsygan 1992; Muslimov & Harding 1997; 
HM98), a simple general expression is not available. For 
the unsaturated case, the gap is pancake-shaped, and the 
boundary conditions at the open field lines are not impor- 
tant. Thus one has (V • E)||/(47r) ~ (V • E)/(47r), which is 
essentially (p — p^j), so that 



/(I) 



Pgj 



rj'^Kil - r? 3) 



(1 - v'i'^v ^) 



3 „-3n 
e,6) 



1 - 0.15r 



(65) 



e,6 



where eqs.(10,ll) of HM98 has been used, (sina) term 
is neglected, and 77* — R/Re = Re q, V = r/RE, n = 



0.15 



{rg/RE){I/MR^) = koRe^s^ «o = {rg/R){I/MR^) 
have been adopted. Since usually S'c(I) <C R, with (|2| 
one can further derive 



■'Note that this estimate of 
been included here. 



/(I) ~ 0.23R^%z, (66) 

is much smaller than that given in eqs. (55,59) of HM98 because the screened accelerator length has 
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where 

Z = S,{\)IRe - 4.2 X 10-2s;4/7p4/7^9/7(^^g^)-3/7^ 

and eq.(^6|) has been used. If this backflow energy is en- 
tirely converted to thermal emission at the hot polar cap, 
then the maximum polar cap thermal X-ray luminosity 
should be 

-^x,pc,max — -^e+,maxj (^^) 

which is 

^ g.lxlO^Vg • s-iB^//2^-9/^i?//'(cosa)5A, 

(69) 

for regime I. The maximum polar cap temperature (as- 

V = 0pcR = 1-45 X 



suming an area of Trr^^, where 
10*P^^/^cm) is hence 



pc,max 



(I) 



(T'Kr'i 



(I) 



1/4 



2.2 X l^fKBl'l^P' 



1/14-2/7 
' e,6 



5.67 X 10 ^ergs • cm 



(cos a 



,5/28 



(70) 
is Stefan's 



where a 
constant. 

For regime II (saturated case), Sc is larger than the ef- 
fective polar cap radius rpc,E = fpcI^E^- The gap shape 
becomes narrow and long. The boundary conditions at 
the open field line boundaries become important, so that 
(V • E)||/(47r) can be much less than (V • E)/(47r). The 
approximation ( |65| ) no longer holds. To achieve a very 
accurate description of such reversed positron fraction, 
one has to appeal to complicated mathematical descrip- 
tions (e.g. Muslimov & Tsygan 1992; IIM98) or numerical 
simulations to get (V • E)||/(47r). However, noticing that 
(V ■ E)||/(47r) ~ (V • E||)/(47r) ~ dE\\/dr, we can stiU get 
an approximate formula for / by adopting (A4) of IIM98, 
and making some reasonable simplication, which finally 
reads 



/(II) ~ 3Ko(i-e' 



-3 



5.7 X 10-VJi?|^eP-i 



(71) 



where the function /(?]) is the correction factor of the po- 
lar cap radius in curved spacetime (see eq.[8] of IIM98), 
which is equivalent to the function q{r) introduced in 
this paper (see eq.|l^]). Other notions have been de- 
scribed above, and typical values have been adopted. With 
eqs.[|8 64|6^j7l[], it is easy to estimate the polar cap heat- 
ing luminosity in regime II, 

ix,pc,max(II) ^ 9.1xl026erg-s-iBp,i2P-i3/V-2i?|g(cosa) 

(72) 

and the polar cap temperature 



5/4 



,(II) ~ 0.53xlO''KBpV2-P 



1/4 p-9/16^-l/2p5/4 



'e,6 



5/16 

(73) 

We see that rpc,max(II) is sensitive to pulsar period (oc 
p-9/16) while Tpc^max(I) is uot (oc P"!/"). This is just 
the refection of the open field line boundary conditions 



on the saturated gap. We also notice that the peak en- 
ergy of such thermal emission falls into the ROSAT band. 
However, it is notable that the thermal X-ray luminosity 
(eqs.[^ or ^^) does not have the dependence that is ob- 
served in ROSAT band {L^{ROSAT) cx L^d, BT97). This 
rules out the possibility of interpreting ROSAT luminosi- 
ties in terms of only thermal emission due to polar cap 
heating. We will show next, however, that a combination 
of such thermal luminosities with the non-thermal lumi- 
nosities due to the ICS of higher generation pairs could 
roughly reproduce such dependence. The thermal com- 
ponent discussed above will dominate the ROSAT band 
emission in most of the millisecond pulsars. Another point 
is that, the luminosity an d p olar cap temperature formu- 
lae presented above (eqs. p9| , [70| , [72|j73| ) are still upper lim- 
its, especially for the saturated regime (regime II), though 
actual values are not much lower. For the millisecond pul- 
sars which are near the deathlines (Zhang et al. 1999), 
such deviations could be even larger, since our approxi- 
mations adopted to describe the accelerator are no longer 
accurate enough. 

3.3. Non-thermal X-ray luminosity 

In our full-cascade picture, we can naturally get a non- 
thermal X-ray emission component from the soft tail of the 
ICS spectra of the higher generation pairs. Attaining an 
accurate spectrum requires a more careful modeling using 
Monte Carlo simulation approach, which will not be done 
in this paper. However, some simple estimate of the non- 
thermal X-ray luminosities of pulsars coul d be ach ieve d 
using the analytic method described in Sect. 2.3 and 2.4. 

Let us compute the luminosity below a certain energy 
Ec- In principle, one should add the contributions from 
all the branches over the whole cascades. However, it is 
notable that the SR spectra of all the SR branches can 
not get down to as low energy as the X-rays observed by 
ROSAT and ASCA (Harding & Daugherty 1999; Rudak & 
Dyks 1999). This could be justified as follows. The lowest 
characteristic SR energy in the full-cascade is 



^c,mi„(SR) = EoKsk'^^'^, 



(74) 



and in principle, the minimum SR energy is the blueshifted 
cyclotron energy (corresponding to the electron occupying 
only one Landau state, see Harding & Daugherty 1999) 



-Bmin(SR) = 7<;sj,j|?iWB,e = 7csr, || mc^B^ , 



where 



7CSR,II 



(CsR-1) 



-7711 



(75) 



(76) 



is the initial parallel energy of the last generation (CsR-th 
generation) pairs. 

For a typical normal pulsar with P = 0.1s, i?p,i2 = 5.8 
■(thus Be,i2 ^ Bp^i2R~e\ = 1 for Rsfi = 1 



a = 30°, we have E^^ - 196GeV (eq. 



and 

_ E^^^ = 0.52GeV 
(eq.[^), ^SR -il/16 (eq.l32|), CsR = 3.14 (eq.Q), and 
r/|| = 0.237 (eq.H), so that i;c,min(SR, nPSR) - 32.8MeV 
and i;min(SR, nPSR) = 1.26MeV. Both values are much 
higher than the X-ray bands we are interested in (e.g. be- 
low lOkeV). Similarly, for a typical millisecond pulsar with 
P - 0.005s, Bpg = 0.5, Re& = 1.2, and a = 30Vwe 
have Eq ~ 163GeV (eq.|l|), E^^^ = 159GeV (eq.pl). 
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KSR - 1/11 (eq-il), CsR = 1.01 (eq.@), and ry,, = 
1.6x10-5 (eq.|§), so that £;c,mi„(SR, msPSR) - 149GeV 
and £;mi„(SR,msPSR) = 0.026KeV. Note although the 
lowest end of this SR spectra could be down to as low as 
0.026KeV, there is a very wide span of energy of the last 
SR generation starting from more than 100 GeV. Thus the 
SR contribution to the ASCA and ROSAT band X-rays is 
tiny compared to the ICS contributions as described be- 
low. 

For the ICS branches, the expression of characteris- 
tic energy (eq.p^) is similar to the minimum SR en- 
ergy (en. |75l), but the physical meaning is different (see 
Sect j2.3| ). The most important difference is, as long as 
the resonant scattering condition (eq.Q]) is satisfied, the 
Lorentz factor of the pairs will keep decreasing so that the 
characteristic energy of the ICS process will cover a broad 
energy range, the low energy end of which could present 
a significant contribution to the X-ray band we are inter- 
ested in. 

As described above, the complication of the full-cascade 
picture makes it difficult to obtain an analytic expres- 
sion which includes all the possible contributions to the 
X-ray luminosity. The main complication arises from the 
"mixed" generations. However, there are two essential fea- 
tures that make such an analytic expression possible: (1) 
The energy distribution to the two branches is fixed for a 
certain pulsar (e.g. ri± and 7/|| are constant for a certain 
pulsar regardless the generation order, see eqs. p7| , ^ ) if 
we ignore the small variation of Be in the different gen- 
erations; (2) Only ICS branches can have significant con- 
tributions to X-ray luminosities. In principle. The upper 
limit of the non-thermal X-ray luminosity below a certain 
energy Ec could be estimated as 



SR-photon-produced pairs. 



int(CsR) 
Lx,nth(Ec) <y Lpc y ^ 
k=l 



int(Cics,k-l) 



(77) 

Here ^sr is just the pure SR generation order parameter 
defined in (Q, and Cics.k is some more general ICS gen- 
eration order parameter defined by 



C: 



ICS,k 



log(Kics) 



+ 1, 



(78) 



which is the number of pure ICS generations for the typical 
energy of the k-th SR generation. 



Ek = EnK 



O^SR' 



(79) 



to reduce to the escaping energy E^sc (eqs.|2^, 29 1). Note 
that we have defined that, for k — 0, Ek is just Eq (eqs.psf 
or 1^). To keep consistency with ( p4[ ) and (p5|), we still 
keep the term (+1) in (78), but usef/c — 1) in ( ]77| ) and 
( |8l| ) to subtract this extra generation. Another important 
parameter in (|77|), ?7c,fcj", is to describe the portion of the 
particles' energy that goes to the band below certain Ec 
in a certain ICS branch, which is defined as the ratio of 
the Lorentz factor required to produce ICS photons with 
energy E^ 

7, = Se/[(l-/?/i)2.8fcr], (80) 

to the initial parallel Lorentz factor of particles from the 
(j — l)-th ICS-photon-produced pairs from the {k — l)-th 



. _ fOKsR_Kics_ 

Ik,] — „ 



SO that 



1, 



Ic > Ik,] 
Ic < Ik,] 



(81) 



(82) 



and the truncation of ?7c,fcj" = 1 at 7^ > 7^^ ensures the 
energy portion can not exceed 100%. Note that ( |77| ) has 
included most of the branches that have contributions to 
the X-ray luminosity. There are still some branches miss- 
ing, e.g., the ICS branch from the SR branch produced by 
the even earlier ICS branch. In any case, the correction of 
these "higher order" branches to ( [77| ) is tiny. 

The estimate in ( [tt] ) actually involves the assumption 
that the conversion efficiency from particle kinetic energy 
to radiation in the ICS branches is 100%. This is only 
true when the resonant scattering condition ^ is ful- 
filled. If the resonant condition fails, the estimate in (77) 
is then just an upper limit. Thus it is essential to com- 
pare 7res (eq. [|) with 7, (eq.|80|). If 7res < 7c, the ICS 
processes producing photons with Ec are well within the 
resonant regime, then rjc,k,] is essentially the precise en- 
ergy portion. If 7res ^ 7c, the ICS processes can still 
produce photons with Ec- But the process just marginally 
occurs in the resonant regime, thus 77c,fe,j is a marginal 
upper limit. If, however, 7103 3> 7c, the resonant scatter- 
ing has stopped at a much higher energy than 7c. Then 
only a small portion of the photons with Ec could be pro- 
duced, then rick.j could just be regarded as a very loose 
upper limit. Such differences are refiected in Table 3. 
Note that when calculating the X-ray luminosity within 
a certain band, Ed < E < Ec2, we have actually taken 



Lai 



ath 



(Ec: 



Lx.nth{Eci) as the predicted lumi- 



nosity in this band. In Table 3, we have adopted 
and '<', for the cases of 7103 < 7ci, 7ci < 7ros < 7c2, and 
7rcs > 7c2, respectively. For old normal pulsars such as 
PSR 1929-f 10 and PSR 0950+08, since both the cooling 
temperature (^) and the polar cap temperature ( [73| ) are 
much lower, the ICS efficiency is lower. Unless 7ros < 7ci, 
we will adopt < since the actual luminosity could be much 
lower than the predicted upper limits. 

The complicated form in equation ( |77| ) could be sim- 
plified for the pulsars in which pair production from ICS 
photons are not important. This is the case for most pul- 
sars except the ones with high B (Table 2). For Cics.k < 2, 
equation (|77|) turns out 



Lx,nth{Ec) < 



int(CsR) 



V\\Vc,k.i)- 



(83) 



fc=i 



In our calculation, we have adopted ( |77| ) uniformly with- 
out losing generality, and used ( p3| ) to do the comparison. 
We found that, for high B pulsars, the luminosity from 
( [TTI ) is much higher than the one from (83), which means 
that the higher order generations have important contri- 
butions to these pulsars. 

The use of (^) for the millisecond pulsars is not ac- 
curate. This is because, for millisecond pulsars, we have 
7c > 71,1 • In other words, the maximum emission energy 
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of the first ICS-branch (i.e. the ICS-branch from the sec- 
ondary pairs) has dropped below the ASCA and ROSAT 
band. This can be understood by noting the very small 
r/||(msPSR) (eq.|39|) and Kics(msPSR) (eq.|4|). Thus us- 
ing ( [77]) will underestimate the non-thermal X-ray lumi- 
nosities. To amend this, we introduce a separate formula 
for millisecond pulsars 

int(CsR) 

L7:Xtf'{Ec)^L,, E ivT'i'^'"')' (84) 

k=l 



where Cc,fc is defined by 



Cc 



0, 

log(-Ec/-Efc 

log KiCS 



Ec < Ek 



(85) 



Ec is the energy we are interested in, and follows ([79|). 
^c,fc is the number of the ICS generations with which the 
typical emission energy is dropped to E^ from Ei^. For 
millisecond pulsars, this value is smaller than unity. Thus 
equation ( |84| ) can raise the X-ray luminosity estimate a 
little bit compared with ( [77| ) or (|3|) . Note again the adop- 
tion of <, < or is pending on the resonant condition. 

Finally, we should notice another general relativistic ef- 
fect, i.e., the redshift of the photons, when comparing the 
theoretical predictions with the observations. More specif- 
ically, the energy of interest E'c, which is used in above cal- 
culations, is connected with the observational energy i?obs 
by a redshift factor: 

i?c = (1 - ^)-'/'i?obs =^ (1 - 0.4r-i)-i/2i?„bs. (86) 

The redshift correction factor (1 — ^)~^^^ is 1.13 for Tg^e = 
1.8 and 1.29 for re,6 = 1- By adopting Kbs = 0.7- lOkeV 
and 0.1-2.4keV, we can then calculate the (maximum) non- 
thermal X-ray luminosities for the ASCA and ROSAT 
bands, respectively. 

3.4. Results 

We have performed detailed calculations of the predicted 
non-thermal and polar cap heating thermal X-ray lumi- 
nosities of all the known X-ray pulsars detected by ROSAT 
or ASCA. In Table 2, we listed the model parameters of 
these 35 pulsars, which are separately grouped into nor- 
mal and millisecond pulsars. Table 3 shows the calculation 
results of both the thermal and non-thermal X-ray lumi- 
nosity predictions of the known X-ray pulsars and the com- 
parison with the observational data (BT97 and S98). The 
model prediction for the total ROSAT band is a sum of the 
ICS-origin non-thermal luminosity (|Q or [ p 4 | ), the ther- 
mal luminosity due to polar cap heating ([|69| or |72]), and 
the cooling thermal luminosity by assuming R — lO^cm for 
all the pulsars. This is because the predicted peak thermal 
components usually fall into the ROSAT band. For non- 
thermal luminosity both for ROSAT and ASCA bands, 
we use (^^ for all the normal pulsars, and use (^) for all 
the millisecond pulsars. Unless the temperature data are 
available from data analysis {Tgfi ~ 1.5 for Vela, Ogelman 
et al. 1993; g ^ 0.56 for Geminga, Halpern & Wang 
1997; T,,6 0.8 for both PSR 0656-1-14 and PSR 1055- 
52, Greiveldinger et al. 1996, Wang et al. 1998), we have 



adopted the values of Tg^ and T^fi according to ( p^) and 
( |70| ) or (^ . For other model parameters, we have adopted 
£, — I (the last open field line), a = 30° for all the pulsars, 
Re.6 — 1-8 for most of the normal pulsars, Rs.e — 1-5 for 
the relatively old normal pulsars (PSR 0950+08 and PSR 
0823-1-26) and Re,6 = 1-2 for the the millisecond pulsars. 
Such adoptions for Re,6, though rough, are based on the 
numerical results of HM98 (their Fig. 10). The weak de- 
pendence of Re,6 on pulsar parameters like P, Bp and T 
is the main reason for us to adopt fixed values oi Re,6 for 
different groups of pulsars. In principle, for normal pul- 
sars, Re,6 is not a free parameter, since it is determined by 
Thfi and the pulsar parameters (P, Bp 12 and cos a), and 
Thfi itself is also the function of P, Bp^i2 and cos a (see 
eqs.0 or |Q) if we assume the polar cap heating is the 
only process to decide polar cap temperature. However, 
the complication of the processes prevents an analytic ex- 
pression of 
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observed value 
predicted upper limit 
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polar cap heating 
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Fig. 4. — Comparison of the luminosities observed and predicted 
in the ROSAT band for 29 pulsars. Note upper limits are adopted 
for some pulsars both for the observational and theoretical values. 
Note how Lx{ROSAT) oc Lg^ feature is well reproduced, and how 
the thermal component due to polar cap heating dominates the lu- 
minosity in millisecond pulsars. 

The X-ray luminosity predictions versus observations 
are plotted in Fig. 4 and Fig. 5. If a pulsed signal is detected 
for a certain source, we have adopted the pulsed emission 
luminosity as the observational value, since we believe that 
only this luminosity is relevant to what we are discussing 
in this paper. For the detections without pulsed signals, 
we adopt their luminosities as upper limits. For the theo- 
retical results, we adopt the sum of both the thermal and 
non-thermal X-ray luminosities for the ROSAT band, but 
only non-thermal luminosity for the ASCA band. We re- 
gard both the values marked with and in Table 2,3 
as the actual values, while the values marked with '<' as 
upper limits to plot Fig. 4 and Fig. 5. It is notable from 
Fig. 4 that, a Lx{ROSAT) oc Lgd feature is roughly re- 
produced, and the thermal component due to polar cap 
heating clearly dominates the ROSAT band luminosities 
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of the millisecond pulsars. The comparison of the theoret- 
ical predictions with the ASCA observations (S98, Saito 
et al. 1998) is less satisfactory, but not in strong conflic- 
tion when one bears in mind the small sample and large 
uncertainties in the ASCA data. Furthermore, the incon- 
sistency could be weakened by adjusting free parameters 
like a and i?£;,6- 
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Fig. 5. — Comparison of the luminosities observed and predicted in 
the ASCA band for 15 pulsars. Notations are same as those adopted 
in Fig.4. 

In Fig. 6 and Fig. 7, we plotted the dependence of our 
model predictions on the free parameters a and Re^- 
We see both the 7-ray and X-ray luminosities are mildly 
decreasing when a gets larger. The decreasing becomes 
prominent when a gets closer to 90°. However, one should 
notice that such decreasing can not get too low (say, be- 
low 2 orders of magnitudes), since the smaller sin a term, 
which has been neglected in writing down the parallel elec- 
tric fields in the HM98 model, will become important for 
a ~ 90°. The dependences of the non-thermal X-ray lumi- 
nosity on REfi or a are not smooth (Fig. 6a, 7a, 7b). This 
is expected since the zigzags actually reflect the contri- 
butions from different ICS branches. For different adop- 
tions of Re,^ of a, the significance of the contributions 
of certain ICS branches to the final luminosity in a cer- 
tain band actually varies. Thus the zigzag features are 
real physical effects. However, we notice that there ex- 
ist several non-differentiable points in the curves. This is 
an artifact, which arises from the truncation features in 
(|8^). We believe that the detailed numerical simulations 
can smooth out these non-differential points, though the 
small "bump" features should remain. The present es- 
timate could bring some errors near these non-differential 
points. For the millisecond pulsars, the height dependence 
is simpler (Fig. 7b), mostly because the generation struc- 
ture is simpler (both ^sR and Cics are smaller than 2), and 
also due to the non- integer description (p4[). 

4. CONCLUSIONS AND DISCUSSIONS 



4.1. Conclusions 

We have modified the conventional CR(or ICS)-SR po- 
lar cap cascade picture by including the contributions 
of the ICS photons produced by the higher generation 
pairs within the framework of the space-charge-limited 
flow acceleration model proposed by Harding & Muslimov 
(HM98). Such a full-cascade picture is more complicated 
than the canonical model. The complications mainly re- 
sult from the two-branch feature of the cascade, and the 
different regimes of the ICS process under different con- 
ditions. The important parameters, i.e., the temperatures 
of both the full surface and the hot polar cap, are still 
uncertain since there are many other factors (besides cool- 
ing and polar cap heating) that can influence them. The 
uncertainty of the temperatures will bring uncertainties of 
the "resonant ICS condition" (eq.[Q), the height of the 
accelerator Re,6, etc., and thus bring uncertainties to our 
model predictions. Furthermore, the complication of the 
accelerator itself (HM98) has already been difficult to de- 
lineate analytically. Therefore, to achieve a more accurate 
prediction, one has to appeal to detailed Monte Carlo sim- 
ulations. However, some interesting features, e.g. the reso- 
nant condition is satisfied for a wide range of parameters; 
Re,6 does not sensitively depend on pulsar parameters; 
the field strength variation is not important within gen- 
erations; only ICS branches have notable contributions to 
X-ray luminosities; and so on, make an analytic descrip- 
tion of the process possible. We have tried to treat the 
process as accurately as possible so as to keep as much 
information as we could. Such efforts include the correc- 
tion of the general relativistic effects (enhancement of the 
field strength and photon energy redshift); the correction 
of the non-zero pitch angle effect for higher SR genera- 
tions; more accurate and separate adoptions of the critical 
X for describing the cascade and escaping processes; sepa- 
rate treatments of the normal and the millisecond pulsars, 
and so on. As shown in the previous sections, our analytic 
description has presented an overview to the high energy 
emission features of the pulsar population. Our main find- 
ings in this paper can be summarized as follows: 

1. When high energy photons are converted to pairs 
in the pulsar magnetosphere via 7 — i? process, not only 
the perpendicular energy of the pairs are further converted 
into radiation via SR, most of the parallel energy of the 
pairs will be also converted to radiation via the ICS of 
the particles with the soft thermal photons. Thus there is 
almost no energy loss during the conversion of the acceler- 
ating particle energy to radiation in the cascade process. 
Such a full-cascade scenario gives a natural interpretation 
to the oc {Lsd)^^^ feature. 

2. SR branches of the cascade have little contribution 
to the X-rays in the ROSAT and ASCA bands. The ICS- 
branch spectra, however, can result in a non-thermal X-ray 
component with their soft tails. Such non-thermal X-rays 
of ICS-origin play an important role to account for the 
ROSAT and ASCA luminosities. 

3. The polar cap heating from backfiowing particles of 
the inner accelerator can contribute a hot thermal com- 
ponent with small area (e.g. polar cap), which also has a 
contribution to the ROSAT luminosity. 

4. The combination of the thermal and the non- 
thermal components in the model reproduces the re- 
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Fig. 6. — The inclination angle (a) dependence of the predicted luminosities, (a) for a normal pulsar: P = 0.1s, 
■Bp, 12 = 5.8, re,6 = 1-8, note ROSAT luminosity is dominated by the non-thermal component; (b) for a millisecond pulsar: 
P = 0.005s, -Bp,i2 = 5 X 10~^, Re,6 = 1-2, note ROSAT luminosity is dominated by the thermal component. 
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Fig. 7. — The emission height {Rsfi) dependence of the predicted luminosities, (a) For a normal pulsar: P = 0.1s, 
■Bp, 12 = 5.8, a = 30°. b) For a millisecond pulsar: P = 0.005s, Bp. 12 = 5 x 10^^, a = 30°. The zigzags indicate 
the contributions from different ICS branches. Note also that the ROSAT luminosity is dominated by the non-thermal 
component for normal pulsars, but by thermal component for millisecond pulsars. 
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ported L^{ROSAT) ~ lO^^L^d feature (BT97). The non- 
thermal components predicted in the ASCA band is not 
in severe contradiction to the ASCA data (S98). 

5. For the ROSAT band, the non-thermal components 
dominate the luminosities for all the normal pulsars, while 
the upper limits of the thermal luminosities are usually 
higher than the non-thermal luminosities for the millisec- 
ond pulsars. We thus predict that the unidentified spectral 
features of the pulsed emission of the millisecond pulsars 
might be of thermal origin. 

Zhang & Harding (1999) have also presented similar 
conclusions. The main improvement of this paper is that 
we have included both the unsaturated and the saturated 
cases of the space-charge-limitcd flow accelerators. 

4.2. Comparisons with the outer gap model 

The Lx ^ lO^-^Lsd feature has been also interpreted 
in terms of the thick outer gap model by Cheng et al. 
(CGZ98; CZ99), thus it is of interest to compare our model 
with the outer gap model. The main differences between 
the two models lie in the following aspects: 

1. Wc interpret the non-thermal X-rays as being due 
to the ICS of upward polar cap cascades, while the outer 
gap model attributes the non-thermal X-rays to the SR of 
the downward cascades from the outer gap particles. A 
prominant feature of our model is that we can simulta- 
neously reproduce both the observed oc (Lgd)^^^ and 
the Lx ~ 10~'^Lsd dependences well. The thick outer gap 
model, however, has a quite different Lgd-dependence of 

the 7-ray luminosities {L^ (x i3p'72-P~^^^> see a combi- 
nation of eqs.[23], [22] of Zhang '& Cheng 1997) which 
severely violates the observed cx (Lgd)^^^ dependece, 
although Lx ~ 10~'^Lsd feature was also well reproduced. 

2. As for the thermal X-rays, our model predicts a hot 
small area component in terms of polar cap heating, and 
adopt the full-surface temperature from a simple cooling 
model. In the outer gap model, both the full-surface and 
the polar cap temperatures are predicted in terms of the 
thick outer gap heating (Zhang & Cheng 1997; CGZ98; 
CZ99). CZ99 also took into account the possible effect of 
the polar cap heating in old pulsars PSR 1929+10 and PSR 
0950-1-08, in which the outer gap might not be formed. An 
important difference is that, for our model, as long as the 
resonant condition (^ is satisfied, our model predictions 
are not sensitively dependent on the surface temperatures, 
while in the outer gap model, the surface temperatures are 
closely related to the non-thermal X-ray luminosities due 
to the self-consistency of the model. 

3. In the pure non-thermal-origin outer gap model 
(CGZ98), some direct analytic expressions of the soft X- 
ray luminosities in terms of the spin-down energy (their 
eqs.[15,16]) are available. However, such a model is in- 
complete since the thermal emission should have a non- 
negligible contribution to the ROSAT band emission, and 
such thermal components have indeed been identified from 
some pulsars (see Sect.l). CZ99 improved the model by 
incorporating the thermal contributions. But they as- 
sumed that the hard polar cap thermal components do 
not have a significant contribution to the ROSAT obser- 
vations. Their results combining both the non-thermal 
and the soft full-surface thermal components can also re- 
produce the Lx{ROSAT) ~ lO^'^Lsd feature. In our polar 



cap model, a simple expression for the hot thermal lumi- 
nosity (p9| or Its]) is available, but no simple analytic ex- 
pression (in terms of pulsar parameters) could be achieved 
to describe the non-thermal luminosity due to complex- 
ity of the cascade processes. A sum of both the thermal 
and non-thermal luminosities can also results in a rough 
Lx[ROSAT) ^ lO^^Lsd dependence. The distinguishing 
point between the two models might be the different ori- 
gin of the soft X-rays from the millisecond pulsars. As 
addressed above, our model predicts that the pulsed emis- 
sion of the ROSAT-observed millisecond pulsars might be 
dominated by the thermal component due to polar cap 
heating, while the outer gap model attributes the millisec- 
ond pulsar soft X-rays also to the non-thermal SR com- 
ponent, although they adopted an additional assumption 
of strong multipole magnetic field components for these 
pulsars. Such an issue could be a discriminator between 
the two models. Unfortunately, among the 10 X-ray mil- 
lisecond pulsars, spectral analyses are only performed on 
two pulsars. Although there is evidence of the existence 
of a hot thermal component in PSR J0437-4715 (Becker 
& Triimper 1993; Halpern, Martin & MarshaU 1996; Za- 
vhn & Pavlov 1998), a combination of the ROSAT and 
ASCA data shows that the spectrum of PSR B1821-24 is 
a power-law, and thus, of non-thermal origin. According 
to our results (Table 3), the observational values both in 
ROSAT and ASCA band for PSR 1821-24 are much higher 
than model predictions, and also much higher than the ob- 
served values of other millisecond pulsars. This leads us to 
suspect that there might be some special non-thermal X- 
ray emission mechanisms operating in this pulsar. Thus, 
when more data are accumulated, the spectral analyses of 
the pulsed emission of the other millisecond pulsars could 
provide evidence for or against our idea. It is worth notic- 
ing that the thermal feature for the millisecond pulsars is 
also expected by Zavlin & Pavlov 1998). However, if there 
indeed exist strong multipole magnetic fields near a mil- 
lisecond pulsar's surface, the differentiating between the 
models becomes difficult, since our model also predicts 
much stronger non-thermal emission components due to 
the enhancement of the ICS emission efficiency. 

4. Among 35 detected spin-powered X-ray pulsars, 
pulsed emission components are only observed in 10 pul- 
sars by ROSAT, 4 of which (Crab, Geminga, PSR B1509- 
58 and PSR B0540-69) are also observed by ASCA. ASCA 
also observed pulsed strong non-thermal emission from an- 
other pulsar, the millisecond pulsar PSR B1821-24. For 
other pulsars, only the total X-ray luminosities are re- 
ported. Our polar cap model actually predicts pulsed 
emission luminosities. Thus, as long as the pulsed luminos- 
ity is detected, we use it to compare with the model pre- 
diction. For those pulsars in which no pulsed emission is 
detected, we regard the total luminosity as an upper limit 
(see Fig. 4, 5), since besides the internal (cooling and heat- 
ing) and magnetospheric (polar cap or outer gap heating, 
and ICS or SR non-thermal emission) origin X-rays, there 
are still many mechanisms, such as the accretion from the 
interstellar medium (Paczynski 1990), pulsar wind nebular 
emission (Arons & Tavani 1993), etc., which can account 
for non-pulsed X-ray emission. We noticed that, however, 
in the outer gap model, the total luminosities are adopted 
to compare with the model predictions. On the theoret- 
ical aspect, the outer gap model predicts a precise value, 
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but our model predictions could either be a actual value 
or an upper limit pending on whether the ICS occurs in 
the resonant regime (see Table 3; Fig. 4, 5). More precise 
predictions will be available after detailed numerical sim- 
ulations. 

5. Besides the spectral behavior of the millisecond pul- 
sars, another issue could also be regarded as a differential 
criterion of the competing models. For the relatively old 
canonical pulsars such as PSR 0950+08 and PSR 1929+10, 
the outer gap model predicts pure thermal emission (due 
to polar cap heating) , since outer gaps could not be formed 
in these pulsars (CZ99). Our model, however, predicts a 
non-thermal plus thermal feature. The thermal emission 
component due to polar cap heating has a luminosity com- 
parable to the ICS-origin non-thermal component, so that 
both of them are detectable. 

4.3. Further remarks 

We have some further remarks on our model. Firstly, 
the analytic description approach adopted in this paper 
has some limit atio ns in presenting detailed predictions. As 
shown in Sect j2.3| , when deducing the recursion relations 
between adjacent generations, we have only examined the 
characteristic emission energies, but ignored the detailed 
radiation spectrum of each generation. This is the main 
reason for the non-smooth features in Fig. 7a. Further- 
more, the expressions of the predicted luminosities may 
need further revisions to include more effects that might be 
important. For example, in (|77|), one may need to incorpo- 
rate a weighting factor to describe the 7-dependence of 7, 
since ( ^n\ ) assumes equal number density for the particles 
with different 7s, which might not be necessarily the case. 
Therefore to get detailed 7-ray and X-ray spectra, care- 
ful numerical simulations are desirable. Nevertheless, we 
can give a rough estimate of the non-thermal X-ray spec- 
tral indices. According to Dermer (1990), the thick-target 
magnetic ICS spectral indices are approximately 1 + r/2 
or r within different energy regimes (see his eqs.[41a,41b]), 
where F is the spectral index of the stable input par- 
ticle spectrum. Monte Carlo simulation results (DH96) 
show that for the higher generation pairs, usually one has 
r ~ 1.5 — 2, so that the final photon spectral indices are 
around 1.75 — 2. This is in good agreement with the ob- 
servations (BT97). 

Our present model is not good at describing the high 
magnetic field pulsars. Although we have tried to ac- 
comendate high i? pulsars by including higher order ICS 
branches (eq.[0), there are still several aspects of high B 
pulsars that we did not take into account. First, the model 
predictions presented in this pulsar are mainly very loose 
upper limits rather than actual values. This is because, for 
high B pulsars, the cyclotron energy is very high. The typ- 
ical Lorentz factor for the resonant ICS condition (eq.[Q) 
is thus also very high (e.g. ~ 10^), and is much higher 
than the Lorentz factor (e.g. about 100) required to emit 
soft X-rays in the ROSAT or ASCA bands. Thus the elec- 
trons may not be able to reduce their energies enough to 
scatter the soft photons to these bands, so that in some 
cases, the pulsar may not be detected in X-rays at all. 
Another point is that in high B pulsars, the accelerator 
should be controlled by ICS, since both the upward and 
downward scatterings occur in the resonant regime, and 
the anisotropy of the ICS processes no longer exist. The 



accelerator is hence moved back to the vicinity of the sur- 
face {Re,6 — 1), and the expression of Eq (Q or 
should be modified by the ICS-controlled one. All these 
will alter the generation structures of the high B pulsars, 
and change the model predictions to some extent. Fur- 
thermore, photon splitting is also important for high B 
pulsars (HBG97), and will complicate the cascade picture 
even more. Thus to compute the cascade spectrum of these 
high B pulsars, one has to appeal to rather complicated 
numerical simulations (see e.g. Baring & Harding 1999). 

Though with limitations, the analytic approach adopted 
in this paper has its advantage to present a general 
overview of the emission properties of different pulsars 
without performing detailed numerical simulations. For 
example, our model shows that, for most young pulsars, 
the non-thermal luminosities are much higher than both 
the full-surface and the hot polar cap thermal luminosi- 
ties, so that the thermal components should be buried by 
the non-thermal power law continuum. For some middle- 
aged pulsars such as Geminga, PSR 0656+14, and PSR 
1055-52, however, the predicted non-thermal X-ray lumi- 
nosities are of the same order of the full-surface thermal 
luminosities, thus these thermal components could be de- 
tectable. This is just what is observed. Another example is 
that our model seems to have the ability to interpret some 
"missing detections". For example, ROSAT has detected 
10 millisecond pulsars. If pulsar X-rays for the millisec- 
ond pulsars are indeed of non-thermal origin as the outer 
gap model argued, then most of them should be detected 
by ASCA as well, since the ASCA's band is higher and 
wider. However, only 2 out of those 10 are detected by 
ASCA. In our model, the ROSAT-band emission for these 
millisecond pulsars is dominated by the thermal emission, 
which is stronger than the non-thermal component. In the 
ASCA band, only the non-thermal component is present, 
so that the predicted luminosities are lower. In fact, our 
predictions of the ASCA luminosities of these 10 millisec- 
ond pulsars are usually too low to be detected (see Table 3) 
except for PSR B1821-24, which has been indeed observed 
to be dominated by the non-thermal emission though our 
prediction is still less by almost an order of magnitude. For 
PSR J0437-4715, we predicted a low luminosity in ASCA 
band, which is close to the observational level. The detec- 
tion of this pulsar by ASCA is due to the close distance of 
the source. Another example of the "missing detections" 
is, among the ASCA detected 16 pulsars (S98), only two 
pulsars, i.e., PSR B1610-50 and PSR 1853+01, were not 
detected by ROSAT. This is not easy to be understood 
in terms of the outer gap model. In our model, we can 
present a simple interpretation. We noticed that these 
two pulsars all have the magnetic fields higher than lO^^G. 
The absence of soft X-rays from these two pulsars might 
just be due to the failure of the resonant scattering con- 
dition in the ROSAT band as discussed above. However, 
some other high B pulsars, e.g. PSR B1509-58 and PSR 
B2334+61, have been detected by ROSAT. Thus numeri- 
cal simulations are desirable to tell the differences in these 
pulsars. 

There is growing evidence that the optical emission of 
the spin-powered pulsars are of non-thermal origin. The 
ICS of the higher generation pairs discussed in this pa- 
per has a spectrum which could extend to as low as the 
optical band. However, the resonant ICS condition fails 
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in this regime, so that we could not present a plausible 
estimate of the optical luminosities without detailed nu- 
merical simulations. Thus whether or not the observed 
non-thermal optical emission is due to such ICS processes 
remains unclear. 

Radio emission from pulsars is commonly believed to be 
of certain coherent origin owing to the very high brightness 
temperatures observed. Thus pair production, which can 
arise various plasma instabilities in the magnetospheres, is 
believed to be the essential condition of pulsar radio emis- 
sion. Our full-polar cap cascade model presents multi- 
component pairs with different energies, which increases 
the possibilities of the "overtaking" two-stream instabil- 
ities that might account for pulsar radio emission. Such 
detailed models deserve constructing. 

4.4. Detailed comparisons with observations 

It is worth presenting a closer comparison of our model 
predictions with the observations of some well-observed 
pulsars. 

Crab and PSR B0540-69: Crab and its twin in the LMC 
are the youngest known pulsars. Both the ROSAT and 
ASCA observations show that they have very high lumi- 
nosity radiation in X-rays, which are 1-2 orders of magni- 
tude higher than our predictions. This leads one to sus- 
pect that there might be some more efficient X-ray emis- 
sion mechanisms operating in these pulsars. However, our 
model could not be ruled out. We noticed that the lumi- 
nosity data reported by BT97 and S98 are the values as- 
suming isotropic emission of the X-rays. Thus the actual 
values should be much smaller if the radiation is strongly 
beamed. For example, these values could be down by an 
order of magnitude if one assumes the radiation is coming 
from 1 steradian solid angle. For Crab and PSR B0540-69, 
the strong non-thermal spectral features make us suspect 
its magnetospheric origin, and thus, the emission could be 
beamed. 

Vela: Vela pulsar is actually weaker in X-rays than 
expected. It was not until 1993 when Ogelman, Finley 
& Zimmermann (1993) first discovered its weak pulsed 
X-ray emission with ROSAT. Its soft X-ray emission is 
dominated by a thermal component with the surface ra- 
dius only 3-4km, which is much smaller than the conven- 
tional value, i.e. 10km. In fact, in BT97 and S98, the 
detected X-ray luminosities are much lower than the ex- 
pected values according to their empirical relationships. 
Recent RXTE (2-30keV) observations (Strickman, Hard- 
ing & deJager 1999) reveal an even lower luminosity in 
this band. Such a feature does not conflict with our model 
prediction, since what we predicted is an upper limit (see 
Table 3). An interesting discovery of Strickman, Harding 
& deJager (1999) is that although the RXTE spectrum of 
Peak 1 of this pulsar joins to the OSSE spectrum smoothly, 
the spectrum of Peak 2 does show the evidence of another 
component, which would extrapolate to within a factor of 
two of the optical spectrum. Detailed numerical simula- 
tions might be able to show whether the extra component 
is of ICS origin. 

PSR B1509-58 and other high B pulsars: We have 
discussed the possible complications in high B pulsars. 
Our present predictions are not inconsistent with the ob 
servations. It is worth noticing that by adopting (83) 
rather than ( |77| ) to do the calculation (thus assuming 



branches do not pair produce), the prediction is nearly 
2 order of magnitudes lower than the observations for 
PSR B1509-58. This indicates that the higher order ICS 
branches actually play an important role in high B pulsars. 

Geminga, PSR B1055-52, and PSR B0656-H14: these 
three middle-aged pulsars are the ones which have de- 
tailed spectral information, and are identified as "cool- 
ing pulsars" , which consist of a cooling thermal compo- 
nent and a non-thermal power law. As shown in Table 
2, our model predictions show comparable luminosities of 
the non-thermal and full surface thermal components for 
these pulsars, and hence can well reproduce the observed 
features. 

PSR 1929+10 and PSR 0950+08: As discussed in 
Sect.l, the spectral nature of these two old pulsars are 
still unclear, though they all probably have a hot polar cap 
thermal component. Since both the full surface tempera- 
ture and the hot polar cap temperature are much lower in 
these pulsars, the ICS conversion efficiency should be con- 
siderably lower. Thus the no n- thermal luminosity could 
be much lower than the upper limit we predicted. As a re- 
sult, most probably, the thermal component from the hot 
polar cap may have comparable luminosity as the non- 
thermal component, so that such hot polar caps could 
be detectable. Though PSR 1929+10 is fitted by a sin- 
gle small area thermal component by Wang & Halpern 
(1997), it could be equally well fitted by a thermal po- 
lar cap plus a non-thermal component spectrum (BT97), 
which is just expected by our model. Future XMM obser- 
vations on these two old pulsars will test our model pre- 
dictions. As discussed in Sect. 4. 2, the thermal plus non- 
thermal feature of these two pulsars might be also adopted 
as a distinguishing criterion for the competing models. 

Millisecond pulsars: as discussed above, our model pre- 
diction for the millisecond pulsars could be a feature to 
distinguish the rivaling models. For PSR J0437-4715 and 
PSR B1821-24, our model seems to have presented a quali- 
tatively correct picture. The fact that a two component fit 
for PSR J0437-4715 requires a small-area thermal compo- 
nent actually supports our idea of the thermal origin of the 
pulsed soft X-rays from millisecond pulsars. We emphasize 
that the thermal feature we have predicted is only for the 
pulsed X-ray emission from millisecond pulsars. Another 
remark is that stronger multipole fields near the surface, 
if any, will make the polar cap model and the outer gap 
model indistinguishable since both models predict the non- 
thermal dominant feature. 

We are grateful to the referee Werner Becker for his in- 
formative instructions and helpful advice, and Alex Mus- 
limov for important discussions and useful comments. 
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APPENDIX 



DERIVATION OF EQUATION (35) 

The point of this appendix was first noted by Cheng & Ruderman (1977), and was adopted in WSL97. However, as 
shown bolow, a shghtly different definition (i.e. A3, A4) makes (A9) an accurate relation rather than an asymptotic form 
for the relativistic cases. 

Suppose the ith generation 7-ray ej produce {i + l)-th generation pair with 7j_|_i = ej/2 and pitch angle ^tB, then 7j_|_i 
can be generally expressed as 

7.+1 = [1 + (p,+i,,|/mc)2 + 2nBf/^, (Al) 

where n is the number of the Landau levels in the direction perpendicular to the field line. For 1, this could also be 
expressed classically as 

7i+i = [1 + (pi+i,|| /mc)2 + (ft+i,^/mc)2]i/2, (A2) 
where pj+i, II =pi+isin6kB and = Pi+icos8is- Define 

7i+i,|| = [l + fe+i,||/mc)2]i/2 (A3) 

and 

7i+i,± = [1 + {pi+i,±/mcf]'/\ (A4) 

we get 

7^+1 =7'+i,|| +7'+i,±-l (A5) 

and 

-v^ - 1 

COS^ 6*1,3 = 2 _ ■ (A6) 

This is a little bit different from Cheng & Ruderman (1977) who defined 7i+i,|| = Pi+i,|| /rrac and 7i+i,_L = Pi+i,±/''nc. 
Note 7i+i,|| is not the parallel energy of the particle we are interested in, since when the particle emits the SR photons, 
it will receive a recoil force so that its momentum is usually changed. 

Now define a Lorentz frame with 7i+i,|| (different from 7j_|_ij|) in which only perpendicular momentum is left. In such a 
co-moving frame, the emission of SR photons will usually bring no recoil to the particle in the parallel direction as long as 
the particle emits photons in the direction perpendicular to the magnetic field, i.e., the particle is in high enough Laudau 
levels. This is true in most cases we discussed in this paper. (Note when the particle is in a very low Landau state, above 
conclusion is not necessarily true, since the SR photons will be emitted in random directions.) Thus such a Lorentz factor 
is identical to the real Lorentz factor of the particle after it emits all its perpendicular energy via SR. A handy Lorentz 
transformation leads to the relation 



1 "2 

1 - 7,;+i 



cos^gkB= , 'i;"' - (A7) 



Compare [A6] with [A7], we finally get 



^lui = ^2 jb ^1 (A8) 



or 



_ 7^+l _ (\Q\ 

^^+^'11 " 7.+!,^ " [1 + (7^1 - 1) sin^ ^kB] V2 ' ^""'^ 



which is the equation (35) in the text. 
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FULL POLAR CAP CASCADE SCENARIO 



Table A1 

Broad-band high energy luminosities (above leV, most luminosity contributions are from the band above lOOKeV, thus 
denoted as L-y) observed (Thompson et al. 1999) and expected. The luminosity for PSR B1046-58 is adopted by that 
above 400 MeV, following Kaspi et al. (1999). Pulsar data are taken from Taylor et al. (1995). 'o' denotes observed value, 
'm' denotes model prediction, 'f denotes the full-cascade model discussed in this paper, and 'c' denotes the canonical 

CR-SR cascade model. 



Pulsars 




P(ms) 


Log B(G) 




L^(o) 


i^(m,f) 


L-y(m,c) 


B0531 + 21 


33.40 


12.88 


38.65 


35.70 


35.54 


34.94 


B0833 - 


45 


89.29 


12.83 


36.84 


34.38 


34.71 


34.22 


B0633 + 


17 


237.09 


12.51 


34.51 


32.98 


33.56 


33.38 


S1706 - 


44 


102.45 


12.79 


36.53 


34.84 


34.57 


34.12 


S1509 - 


58 


150.23 


13.49 


37.25 


35.21 


34.85 


32.33 


B1951 + 


32 


39.53 


11.99 


36.57 


34.40 


34.64 


34.58 


B1055 - 


52 


197.10 


12.34 


34.48 


33.79 


33.56 


33.45 


B1046 - 


58 


123.65 


12.84 


36.30 


34.34 


34.46 


33.96 



Table A2 



Model parameters of the spin-powered X-ray pulsars. Pulsar data are taken from Taylor et al. (1995). The 4th column 
'regime' denotes the regime; of the accelerator: T indicates the unsaturated case, while 'IF indicates the saturated case. 
Normal pulsars and millisecond pulsars are grouped separately. 



Pulsars 


P(ms) 


logB(G) 


regime 


/ 












CsR 


Cics 


V± 


V\\ 


B0531 + 21 


33.40 


12.88 


1 


0.00010 


6.25 X 


10- 


-'2 


9.66 X 


IQ- 


-3 


3.44 


2.46 


0.669 


0.331 


B0833 - 45 


89.29 


12.83 


I 


0.00019 


6.25 X 


10- 


-2 


7.64 X 


10- 


-3 


3.21 


2.26 


0.704 


0.296 


J1811 - 1926 


64.66 


13.11 


I 


0.00011 


6.25 X 


10- 


-2 


2.58 X 


10- 


-2 


3.41 


2.83 


0.485 


0.515 


J1617- 5055 


69.33 


12.80 


1 


0.00017 


6.25 X 


10- 


-2 


6.71 X 


10- 


-3 


3.25 


2.25 


0.722 


0.278 


B0633 + 17 


237.09 


12.51 


II 


0.00042 


6.25 X 


10- 


-2 


1.74 X 


10- 


-3 


2.62 


1.71 


0.856 


0.144 


B1706 - 44 


102.45 


12.79 


I 


0.00021 


6.25 X 


10- 


-2 


6.56 X 


10- 


-3 


3.16 


2.19 


0.725 


0.275 


B1509 - 58 


150.66 


13.49 


I 


0.00011 


6.25 X 


10- 


-2 


9.68 X 


10- 


-2 


3.42 


3.87 


0.182 


0.818 


B1951 + 32 


39.53 


11.99 


1 


0.00035 


6.25 X 


10- 


-2 


1.62 X 


10- 


-4 


2.97 


1.63 


0.956 


0.044 


B1046 - 58 


123.65 


12.84 


1 


0.00022 


6.25 X 


10- 


-2 


8.07 X 


10- 


3 


3.15 


2.24 


0.696 


0.304 


B1259 - 63 


47.76 


11.82 


I 


0.00049 


6.25 X 


10- 


-2 


7.43 X 


10- 


-5 


2.85 


1.54 


0.970 


0.030 


J0537 - 6906 


16.11 


12.26 


I 


0.00015 


6.25 X 


10- 


-2 


5.94 X 


10- 


-4 


3.30 


1.86 


0.916 


0.084 


B1823 - 13 


101.45 


12.75 


I 


0.00023 


6.25 X 


10- 


-2 


5.26 X 


10- 


3 


3.14 


2.13 


0.753 


0.247 


B1800 - 21 


133.61 


12.93 


I 


0.00021 


6.25 X 


10- 


-2 


1.20 X 


10- 


-2 


3.18 


2.36 


0.635 


0.365 


B1929 -1- 10 


226.52 


12.02 


II 


0.00042 


6.25 X 


10- 


-2 


1.50 X 


10- 


-4 


2.21 


1.38 


0.958 


0.042 


B0656 + 14 


384.89 


12.97 


II 


0.00026 


6.25 X 


10- 


-2 


1.34 X 


10- 


-2 


2.70 


2.09 


0.615 


0.385 


B0540 - 69 


50.37 


13.00 


I 


0.00011 


6.25 X 


10- 


-2 


1.59 X 


10- 


-2 


3.41 


2.61 


0.583 


0.417 


J1105-6107 


63.19 


12.31 


I 


0.00031 


6.25 X 


10- 


-2 


7.01 X 


10- 


-4 


3.03 


1.77 


0.908 


0.092 


B0950 + 08 


253.07 


11.69 


II 


0.00024 


6.25 X 


10- 


-2 


4.39 X 


10- 


-5 


1.55 


1.15 


0.977 


0.023 


B1610-50 


231.60 


13.33 


1 


0.00017 


6.25 X 


10- 


-2 


5.83 X 


10- 


-2 


3.25 


3.20 


0.295 


0.705 


J0538-h2817 


143.16 


12.16 


11 


0.00070 


6.25 X 


10- 


-2 


3.56 X 


10- 


-4 


2.65 


1.58 


0.935 


0.065 


B1055 - 52 


197.10 


12.34 


II 


0.00050 


6.25 X 


10- 


-2 


7.66 X 


10- 


-4 


2.59 


1.61 


0.904 


0.096 


B0355 -1- 54 


156.38 


12.22 


II 


0.00064 


6.25 X 


10- 


-2 


4.65 X 


10- 


-4 


2.65 


1.59 


0.925 


0.075 


B2334 + 61 


495.28 


13.29 


11 


0.00020 


6.25 X 


10- 


-2 


4.96 X 


10- 


-2 


2.81 


2.67 


0.334 


0.666 


B0823 + 26 


530.66 


12.29 


II 


0.00012 


6.25 X 


10- 


-2 


7.41 X 


10- 


-4 


1.59 


1.23 


0.906 


0.094 


B1853-I-01 


267.40 


13.18 


II 


0.00038 


6.25 X 


10- 


-2 


3.30 X 


10- 


-2 


3.11 


2.71 


0.430 


0.570 


J0437 - 4715 


5.75 


9.06 


11 


0.0081 


9.38 X 


10- 


-2 


5.74 X 


10- 


-lU 


1.49 


1.05 


0.9999 


0.0001 


B1937-F21 


1.55 


8.90 


II 


0.0415 


9.38 X 


10- 


-2 


7.22 X 


10- 


-10 


2.50 


1.17 


0.9999 


0.0001 


B1821 - 24 


3.05 


9.65 


II 


0.0216 


9.38 X 


10- 


-2 


2.46 X 


10- 


-8 


2.74 


1.23 


0.9996 


0.0004 


J2124- 3358 


4.93 


8.67 


II 


0.0070 


9.38 X 


10- 


-2 


3.87 X 


10- 


-11 


1.05 


1.01 


1.0000 


0.0000 


B1957 + 20 


1.60 


8.52 


11 


0.0367 


9.38 X 


10- 


-2 


9.66 X 


10- 


-11 


2.05 


1.11 


1.0000 


0.0000 


J1024-0719 


5.18 


8.79 


II 


0.0074 


9.38 X 


10- 


-2 


9.26 X 


10- 


-11 


1.19 


1.02 


1.0000 


0.0000 


J1744- 1134 


4.07 


8.58 


II 


0.0091 


9.38 X 


10- 


-2 


3.13 X 


10- 


-11 


1.14 


1.01 


1.0000 


0.0000 


J1012-F5307 


5.25 


8.74 


II 


0.0068 


9.38 X 


10- 


-2 


5.83 X 


10- 


-11 


1.09 


1.01 


1.0000 


0.0000 


J0218 -1-4232 


2.32 


8.94 


II 


0.0263 


9.38 X 


10- 


-2 


7.42 X 


10- 


-10 


2.21 


1.14 


0.9999 


0.0001 


J0751 + 1807 


3.47 


8.53 


II 


0.0116 


9.38 X 


10- 


-2 


3.17 X 


10- 


-11 


1.26 


1.03 


1.0000 


0.0000 
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Table A3 

X-ray luminosities observed and expected in the ROSAT band and ASCA bands. ROSAT data follow BT97, and ASCA 
data follow S98. The middle 6 columns are for the ROSAT band, and the right 3 columns are for the ASCA band, 'o' 
denotes observed value and 'm' denotes model prediction, the superscript 'tot' indicates total luminosity detected, 'pul' 
indicates the luminosity of the pulsed portion, 'nth' denotes non-thermal, and 'thl' and 'th2' denote thermal emission of 
full surface and polar cap, respectively. 



Pulsars 



L 



sd 



ROSAT 



band 



' nth 



ASCA 



band 



(m) ii°'(m) 



ix(m) 



B0531 + 21 
B0833 - 45 
J1811 - 1926 
J1617- 5055 
B0633 + 17 
B1706 - 44 
B1509 - 58 
B1951 + 32 
B1046 - 58 
B1259 - 63 
J0537- 6906 
B1823 - 13 
B1800 - 21 
B1929+ 10 
B0656+ 14 
B0540 - 69 
J1105-6107 
B0950 + 08 
B1610 - 50 
J0538 + 2817 
B1055 - 52 
B0355 + 54 
B2334 + 61 
B0823 + 26 
B1853 + 01 



38.65 
36.84 
37.97 
37.22 
34.51 
36.53 
37.25 
36.57 
36.30 
35.92 
38.68 
36.45 
36.35 
33.60 
34.58 
38.17 
36.39 
32.75 
36.19 
34.69 
34.48 
34.66 
34.80 
32.66 
35.63 



35.98 
32.70 



31.10 
33.15 
34.29 
33.44 
< 32.11 
32.95 

33.39 

33.06 
30.00 
32.98 
36.21 

29.35 

32.74 
33.42 
31.96 
31.86 
29.83 



35.85 
31.7 



30.62 



34.10 



29.5 

32.15 

36.1 



32.57 



32.49 
32.96 
32.45 
32.17 
31.24 
32.03 
32.45 
31.72 
32.01 
31.00 
32.25 
32.00 
32.05 
29.07 
31.86 
32.44 
31.81 
27.56 
32.18 
30.46 
31.86 
30.54 
31.89 
28.67 
32.01 



J0437 - 4715 

B1937 + 21 
B1821 - 24 
J2124 - 3358 
B1957 + 20 
J1024-0719 
J1744 - 1134 
J1012 + 5307 
J0218 + 4232 
J0751 + 1807 



^ 31.55 
^ 30.98 
<> 31.24 
^31.11 
<> 30.18 
^ 30.89 

< 30.88 
^31.19 
^ 30.80 
^ 31.04 
^ 31.78 
^ 30.88 
<, 30.78 

< 29.72 

< 29.95 
^ 31.35 
^ 31.02 
<; 28.88 
^ 30.59 
^ 30.55 
<; 30.26 
^ 30.48 
^ 29.93 
^ 28.45 
2 30.69 



< 34.25 

< 33.35 

< 34.23 

< 33.73 

< 31.72 

< 33.46 

< 33.74 
^32.18 

< 33.34 
<, 31.92 

< 34.24 

< 33.45 

< 33.31 

< 31.63 

< 32.24 

< 33.98 

< 33.34 

< 30.28 

< 33.59 

< 31.81 

< 31.49 

< 31.82 

< 32.15 

< 30.74 

< 32.83 



< 34.26 

< 33.50 

< 34.24 

< 33.74 

< 31.86 

< 33.48 

< 33.76 
^ 32.34 

< 33.36 
<, 32.02 

< 34.25 

< 33.46 

< 33.33 

< 31.64 

< 32.39 

< 33.99 

< 33.35 

< 30.30 

< 33.61 

< 31.85 

< 32.03 

< 31.86 

< 32.34 

< 30.75 

< 32.89 



37.02 
33.28 



29.79 
32.83 
34.60 
33.79 
32.96 



35.83 
< 31.83 



29.56 

< 32.73 
34.80 

< 32.93 

< 32.78 



30.08 

36.92 

29.00 
34.36 



< 29.92 

< 31.23 
36.11 

< 29.00 

< 34.36 



30.45 < 30.43 



33.15 < 33.15 



< 34.83 

< 33.96 

< 34.63 

< 34.21 
^ 32.33 

< 33.91 

< 34.34 
^ 32.78 

< 33.80 
<, 32.53 
^ 34.45 
^ 33.83 

< 33.80 

< 31.70 

< 32.83 

< 34.58 
^ 33.41 

< 30.84 

< 34.12 
^ 32.42 
^ 32.10 
<. 32.43 

< 32.75 
^ 31.32 

< 33.33 



34.07 
36.03 
36.35 
33.55 
35.20 
33.71 
33.70 
33.58 
35.40 
33.88 



30.86 
< 32.10 
33.24 
30.35 
31.93 
29.48 
29.30 
30.20 
32.75 
31.60 



30.3 



29.8 



23.64 
25.27 
27.09 
22.33 
23.70 
22.73 
22.30 
22.50 
24.72 
22.37 



^ 31.30 
^ 33.13 
2 32.93 
^ 30.99 
^ 32.66 
<, 31.09 
^ 31.20 
^ 30.99 
<, 32.57 
^ 31.41 



29.82 
31.22 
31.54 
29.49 
30.78 
29.59 
29.54 
29.52 
30.96 
29.60 



^31.31 
<, 33.13 
^ 32.95 
^ 31.01 
^ 32.67 
^31.10 
^31.21 
<, 31.00 
^ 32.58 
^ 31.42 



30.08 
33.81 



< 30.08 



32.97 



30.06 
31.46 
31.77 
29.74 
31.03 
29.83 
29.78 
29.77 
31.20 
29.85 



